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CHAPTER I INTRODUCnCN 


1.1 Background 

Mlcrowavo paxamatrlc amplifiers, the subject of this investigation, 
are a class of amplifiers which utilize either nonlinear reactances, or 
reactances which can be varied as a function of time by applying a suit- 
able pump source. The time variation of a reactive parameter can create 
the equivalent of a negative resistance in a certain frequency range, 
and this negative resistance can be used to provide amplification. This 
Is the origin of the term negative resistance parametric amplifier, or 
simply parametric amplifier. There are two distinctive features of para- 
metric amplifiers that are worth mentioning. First, electronic amplifiers 
'nse energy from an electrical source to Incn^ase the power In a deslrtd 
signal waveform, for most amplifiers, such as vacuum tube or transist.or 
amplifiers, this electrical source Is a direct-current (dc) source, while 
for parametric amplifiers. It Is an alternating-current (ac) source. The 
second feature of a parametric amplifier is its capability of low noise 
amplification. A vacuum tube or a transistor Is essentially a nonlinear 
resistor, and It Is well known that any resistor at non-zero temperature 
will generate thermal noise . On the other hand, a parametric cusplif ler 
utilizes mainly a nonlinear reactance, aind a reactance does not contrib- 
ute thermal noise. This was first pointed out by van der Zlel In 19^8 
when he analyzed the mixing properties of nonlinear capacitances [^l2. 

It is this feature that makes parametric amplifiers the moat attractive 


- 1 - 


ORlvlINAl. PAGE IS 
OF ?OOW 


candidate for low noise front ends in coemunication systems. VThen oper- 


ated cryogenically , the noise figure of a parametric amplifier is compa- 


rable to that of a maser, yet its bandwidth and stability are far supe- 


rior to those of a maser. 


Although van der Zlel was the first person to point out the potential 


use of nonlinear capacitance as low noise amplifier, parametric amplifi- 


cation was t.oeoretically shown to bo possible by Faraday in as early as 


1831 [2], and later by Lord Raleigh in I863 [33* I t was, however, one 


hundred years later, when the parametric effect was experimentally observed 


In an electro-mechanical nonlinear capacitance device [^j. In 1957', the 


first realization of a microwave parametric amplifier was finally made 


by Weiss [i], following the earlier prop''?.**'! by Suhl [6], suggesting the 


’ise of the nonlinear effect in ferrites. This caiised widespread Interest 


aiaong microwave engineers, and in the following few years, with high q,’Jal- 


ity semiconductor junction diodes (often referred as varactor diodes, or 


simply varactors) more readily available , semiconductor parametric ampli- 


riers were soon developed through the efforts of many researchers. 


The semiconductor junction diode has a nonlinear caoacltance . If a 


pump source at frequency f^ and a f^nall amplitude signal at frequency f^ 


are applied simultaneously, the nonlinear capacitance behaves like a time- 


varying linear capacitance at f_. The mixing of f. and f will generate 

S ]p s 


a third frequency component, f -f . This frequency Is usually called the 

? 3 


idler freq’uency, f^. The idler frequency is an inevitable by-product cf 


r 


parametric amplification, suppressing it would also suppress the desired 




ajnpllficatioa at f . It shotild be pointed out that the closer the signal 
frequency Is to half the pump frequency, the closer the Idler frequency 
is to the signal frequency, and the more difficult it is to separate 
signal and idler frequencies by filtering. If the signal and idler fre- 
quencies are so far apart that the signal circuit does not pass the idler 
frequency, the amplifier is called a nondegenerate amplifier. On the 
other hand, if the signal and idler frequencies aire very close or if their 
spectra overlap, the signal circuit can no longer distinguish between 
them, and the amplifier is then called a degenerate amplifier. 

• For degenerate auapliflers, the ordinary concepts of noise figure do 
not apply. Degenerate amplifiers are not amplifiers in the usual sense, 
because they gj,/e output at frequencies not Included in the Input. While 
the noise figure of a nondegenerate amplifier is uniquely defined, the 
noise performance of a degenerate amplifiev depends on the type of input 
signal, the type of detector used, and the interpretation of the detector 
output. Still, two kinds of noise figure, single - s.'’ de band and double- 
side 'cand, are often quoted by manufacturers. 

The single-sideband noise figure is used in operation where the input 
signal spectrum Is confined to one side of half-pump frequency. Although 
t.he signal circuit can not distinguish between the signal and the idler, 
a sharp filter in a subsequent stage can be used to select This type 

of operation is characterized by the reduction of 'useful amplifier band- 
width and a certain def;ree of degralition in signal to noise ratio. In 
operation in which the input s'pectr'Ji surro’unds the half-pump frequency, 


the double- side bend noise figure is used. For amplifiers built with 
diodes of the same quality, the noise figure of a nondegenerate amplifier 
is higher than the double- side band noise figure, yet lower than the 
single- si da band noise figure of a degenerate amplifier. However, it 
must be kept in mind that a degenerate amplifier and a nondegenerate 
amplifier can not be comx>ared by their respective noise figures unless 
the system into which the amplifiers are to be incorporated is first 
specified.- 

It should be obvious that in some instances it is possible to real- 
ize system sensitivities calculated from the double -side band noise figinre 
Vhen this is so, the degenerate amplifier would be no doubt the better 
choice. Sven in applications in which single- sideband ncise figure must 
be used, there may well be practical considerations which would stake the 
degenerate amplifier a better choice. By eliminating the idler circuit, 
the degenerate amplifier is a much simpler device to build. Its puap 
frequency is relatively low as compared to that of a nondegenerate ampli- 
fier. Also, as a consequence of circuit simplicity, the broadbanding of 
a degenerate amplifier is easier. 

1.2 Statement of Problems 

Since its initeption, the parametric amplifier has been plagued by 
the problem of having very narrow bandwidth. Numerous researchers have 
proposed solutions for this problem L7ZE9lC9l[lO^Lll-[i2Z. Iii aost cases 
over simplified a-ssimiptions were made and parasitic elements of tne diode 



togetl;ei with signal circuit loss wero neglectsd. This caused significant 
discrepancies in theoretical and actual responses. 

Seidel and Hermann C?! appear to be the first to attempt broadbanding 
the parametric amplifier by use of a multiple-resonator matching circuit. 
They gave design criteria for a filter circuit of a degenerate amplifier, 
using the approach of setting the derivatives of the gain function equal 
to zero at midband. However, the varactor model Is too simple, and loss 
in signal circuit Is not included. 

Matthaei L3J subsequently derived the gain expressions sul cable for 
vrldeband design, using a complete '/aractor equivalent circuit, and demon- 
strated that, by using proper f 11 cars In signal and idler circuits, frac- 
tional bandwidth of 10^ (at r gain of 15 dS) can be obtained. However, 
no direct way of choosing thc^ proper filters is given, and a considerable 
amount of experimenting is involved. 

Kuh and ?ukada [9j developed an approximate synthesis technique besed 
on more rigorous network concepts. Starting from 3ode's theorem on reflec- 
tion coefficient limitation, equations for gain- bandwidth product are 
derived. ?rom a designer's point of view, Kuh and Fukada's tectuiique is 
more tractable than that of Matthaei, but it also suffers from several 
deficiencies. Diode parasltlcs are not Included in the eqxii valent circuit 
and the circulator is assumed to pass both the signal atnd the idler fre- 
quencies. lu lIO]] subsequently derl*/ed a more exacting synthesis technique 
which, altho’Jgn it is more flexible and more elegant, is less tractable 


than chat of Kuh and F'ukada, and also suffers from the same deficiencies. 
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Perhaps the most widely used synthesis technique is that developed 
by DeJager [Hj, with extensions by Connors [l33 Porra and Somerv-io 
L14'3. However, this technique contains numerous approximations , employs 
a rather simplified varactor equivalent cixcuit, and is limited to only 
double tuned si gnal circuits. Sgami [l2] later developed a new design 
theory based on slope parameter concepts. He has included diode para- 
sltlcs and signal circuit loss in his derivations. While the technique 
is more exact, it is also less tractable and is again limited to double 
tuned signal circuits only. 

Many, if not all, of the above mentioned deficiencies can be removed 
if computer-aided design techniques are fully utilized in parametric am- 
plifier design. With the high speed capability of a digital computer, 
one can afford to use i.’ more realistic varacror equivalent circuit, to 
Include signal circuit loss In the computation of amplifier performances, 
and to explore more complicated topologies. In this research report, 
a computer program will be developed for parametric amplifier design 
with special emphasis on a degenerate parametric amplifier in the fora 
of a microwave integrated circuit. 

Because of their ad'/antages of light weight, low cost and mass 
producibllity, microwave integrated circuits have been Increasingly used 
to replace coaxial cables and waveguides in many microwave instruments. 
Along with these advantages, microwave integrated clorcuits also bring to 
microwave engineers some very cnallenging design problems. Among these 
the most serious is perr-aps the lack of tunability . To overcome this 


difficulty, aiicroHave integrated circuits are some tines first built auid 
tested with adhesive copper foil on oversized substrate at relatively 
low frequencies, usually a few hundred negahertz (MHz) . Tuning is pos- 
sible, though very cumbersome, on this low frequency model. Once the 

circuit is tuned to achieve the desired performance, frequency scaling 

* 

techniqiies are used to bring the operating frequency into the microwave 
or millimeter-wave frequency range, usually from a few gigahertz to tens 
of gigahertz (GHz) . This technique proves to be very useful for passive 
network design _15^« However, for active networks, the scaling is far 
more difficult and less satisfactory because of the difficulty in obtaining 
two active devices ha-zing all their parameters (size, junction and package 
capacitances, bulk resistance, and lead inductance) related tc sach other 
by the same factor. This is particularly true for active networks with 
low stability and high sensiti'zity , such as parametric or tunnel diode 
amplifiers. This may well be the reason that only very few microwave 
integrated circuit parametric amplifiers are reported In the literature. 

1.3 Objectives and Outline cf the Present Study 

The objecti’z-e of this study is to da-zelop a computer program for 
parametric amplifier design with special emphasis on practical design 


♦ All circuit dimensions are reduced by a certain factor. Hence, the 
operating frequency is also Increased by the same factor. For more 
details, see l 15 _. 
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consldarations for aiicroHave Integrated circuit degenerate amplifiers. 

To attain this objective, precision measurement techniq.ues must be 
developed to obtain a more realistic varactor equi'/alent circuit, existing 
theory of parametric amplifier must be mc-tified to inclxide the uew equi- 
valent circuit, and microwave integrated circuit properties, such as 
loss characteristics and circuit discontinuities, must be investigated 
thoroughly . 

In Chapter II the basic theory of semiconductor PN junction is 
briefly reviewed. Lumped-element equi’ralent circuits of packaged varactor 
for various frequency ranges are then proposed following a close exami- 
nation of nhe structure of a "ypical varacror package. Techniques for 
precision measurement of driving-point impedances are given together 
with methods for extracting varactor parameters from measured impedance 
data. 

Chapter HI is devoted to the analysis of parametric amplifier cir- 
cuits. The behavior of a varactor under the Influence of a pump source 
is investigated. This is followed by the formulation of nhe small- 
signal immittance matrix of a pmped ’/aractor. Gain and noise flgxre 
expressions for amplifier circuits employed a complete equivalent circuit 
are presented in a manner which makes them suitable for implementation 
in computer-aided design program. 

The study of microwave integrated circuit properties is covered in 
Chapter Tf . Computational methods for characteristic impedance, effective 
dielectric constant, and attenuation constant are given along with n’Laner- 
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leal results for several frequently used substrates. Glrcvilt disconti- 
nuities, such as open circuits and T- jxmetions, axe discussed in detail. 
Analysis and synthesis methods for one particular circuit component, the 
parallel-coupled band-pass filter, are presented. 

The models and calculation methods developed in Chapter II through 
Chapter IV are used in Chapter V to design and construct a 5-5 GHz 
degenerate amplifier. The computer program used for this design is 
described in detail. Power gain and noise figure of this amplifier are 
reported. 

In Chapter VI the results of this study are summarized and sugges- 
tions are given for further research into all aspects of this study. 


CHAPTER II CHARACTSRIZATION OF MICROWAYE VARACTOR DIODES 


2.1 I ntr oduct Ion 

The varactor diode is a semiconductor p-n junction which is generally 
used not for its rectifying properties hut rather for its voltage depend- 
ant nonlinear capacitance provided primarily hy the depletion layer of the 
junction. 3y specifying the impurity profile throughout the junction 
region, the dependence of the depletion width and hence the nonlinear 
depletion capacitance on applied voltage can be controlled to suit the 
intended application. As depicted in Fig. 2.1. (a), the p-n junction is 
formed by diffusing p-type impurity atoms (e. g. boron) into an n-type 
epitarial layer which is grown on top of an n"^ substrate, '^e substrate 
is purely for mechanical support and is heavily doped ( usual I y -iith arsen- 
ic) to reduce the undesired bulk resistance. An ohmic contact is made to 
a anaH circular area on top of the wafer, and most of the epitaxial layer 
is etched away, except that which is directly underneath the contact. 

In this way a mesa of desired geometry can be formed. The substrate is 
then bonded electrically and mechanically to a mounting post, or pedestal 
which is raised from one of the two conducting end caps of the package 
as shown in Fig. 2.1(b). The top of the wafer is connected to the other 
end cap 'using one or more lead wires or bonding straps. The two end caps 
are brazed to a ceramic casing to ensure hermetic sealing. While such 
packages are rugged and convenient for handling, they are atlso a source 
of parasitic elements which become significant at microwave freq’uencies. 
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In this chapter, the electrical characteristics of a vaxav^or will 
be dealt with first. 3y closely eramining the varactor properties and 
the package structxire, a realistic equivalent circuit is then proposed. 
Finally, means for detejalning the equivalent circuit are discussed. 

This includes the meaiaiirement setup, measurement technique, computational 
methods, and a computer-aided optimization technique for obtaining numer- 
ical values for the eqxiivalent circuit elements from the experimental 
data. 




2.2 Circuit Slements of Packaged Varactor Diodes 

2,2.1 Junction Capacitance . Although the formation of a p-n j-;c- 
tion is actuallj done by a diffxision process, let us visualize vfhat would 
happen if two regions of semiconductor material possessing different type 
of conductivity, one of p-type and the other of n-type, are brought into 
contact. Because of the concentration gradient, electrons would diffuse 
from the n-type region into the p-type region and quickly recombine with 
holes. Similarly, holes would diffuse from the p-type region into the 
n-type region and recombine with electrons. This process would leave a 
net positive charge In the previously neutral n-reglon near the junction, 
and a net negative charge in the p-region due to the ionized donor and 
acceptor atoms respectively. As the diffusion proceeds, an electric 
field is set up which retards and finally stops the diffusion of majority 
carriers across the junction. After equilibrium is established, a narrow 
region called the depletion layer or space-charge layer is left ar the 


k 
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j\inctlon, vfhich is swept completely fi-ee of charge carriers by the elec- 
tric field. The potential difference across the depletion layer is 
called the contact potential or barrier potential. 

The width of the depletion layer depends on both the type of semi- 
conductor material and the impurity distribution near the junction. 
Knowing the impurity distribution, the depletion layer width cau be 
found by solving the one dimensional Poisson's equation for the scalar 
potential 'f'(x) 

_ 

with appropriate boundary conditions. Figure 2.2 shows the space-charge 
distributions of the two most commonly treated junction, namely, the 

abrupt junction and the linearly-graded junction. For the abrupt 

♦ 

junction, the width of the depletion layer is given by 


LliEl 


( 2 . 1 ) 


_ 2€(<t -f Y) 

q 



( 2 . 2 ) 


where 

q = the electronic charge (].6 x 10”^^ C), 
<p= the contact potential (volts), 




* The deri’/atlon of Sq. 2.2 is well cover***^ in essentially every semi- 
conductor device theory bock. See, for example. Chapter 5 of 


P(x) = 


qNi 

(a) 



(b) 


Space-Charge Distributions 

(al abrupt Junction 

(b) linearly-graded Junction 


Fig. 2.2 


€ - the permittivity of the semiconductor material (F/m), 

7 = the applied voltage (volts), 

N = the acceptor concentration (m 

3L 

= the donor concentration (m . 

Since varactors are usually operated in the reverse hias region, for 
convenience, a sign convention is adapted such that a reverse bias voltage 
is positive, while a forward bias voltage is negative. Equation 2.2 can 
be written as 

W = W^(l + 7/^)^/" (2.3) 

where 7^ represents the width of the depletion layer at zero bias. The 
junction capacitance is 


G 


JV 



C,, 


(1 + 7/<p)^/2 


(2.4) 


where A is the junction area, and the capacitance at 

zero bias. 

Similarly, for linearly- graded junction, the width of the depletion 
layer at zero bias is 


W 

0 


/ 12 € 4 > 


{ 2 . 5 ) 


where g is the impurity gradient as shown in Fig. 2.2(b). When bias 
voltage 7 is applied to the junction, the width of the depletion layer 


is changed to 
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W = W^(l + 

The junction capacitance Is then given as 


G 


Sf 



jO 

(1 + 7/<p)^^^ 


( 2 . 6 ) 


(2.7) 


In practice, the impurity distribution is far more complicated and 
usually can be appro^dLmated by either a gaussian or a complimentary 
error function, depending on the diffusion process. In this case, the 
capacitance- voltage relationship must be calculated numerically [l7j. 
However, as a first order approximation, the depletion layer width and 
the junction capacitance of a reaLl p-n junction can be expressed 

V = W^(l + V/<^>)^/° (2.8) 



‘J0_ 


I7H 


(1 + V/<f>) 


(2.-9) 


where t.he val'je of n is in the range of 2 to 3* 

2.2.2 Series Resistance . The series resistance of a packaged diode 
consists of two partsi fixed resistance, 3^, and variable resistance, 3^. 
The •^ixed resistance includes that of the two end caps, the mounting post, 
the substrate, and the bonding straps, and is independent of the bias 
voltage. 3^ can be evaluated if the resistivities of the substrate and 
the conductors are given together with the package configuration. The 
•/arlahle resistance is due to the bulk resistance of the semiconductor 
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nesa excluding the depletion layer. For the case of .in epitaxial di')Ua, 
this resistance is dominant and is gi'/en by 

Pp(x) dx + F^(x) dx ] (2.10) 

whore the integration limits for the first term are from the top of the 
epitaxial layer to the ei'ge of the depletion region in the p- region, and 
from the epitaxial layer - substrate interface to the edge of the depletion 
region in the n-region for ti.e second term. The p-type and a- typo resis- 
tivltiBs are functions of the acceptor and donor Impurity concentrations, 
respectively. is thus a function of depletion width which in turn is 
a function of the bias voltage. Consequently, the total series resl tamce 
R is a function of bias voltage. Larger bias voltage causes a wider 

3 

depletion layer, which lowers the series resistance. 

For abrupt junction, 3q. 2.10 can be calculated by niaki.ng seme sim- 
plifying approximations. For instance, if the p-region is assumed to be 
negligibly thick, and f « ? , so that V ^ W , then Sq. 2.10 becomes 

w w • p Q 

where t l.s the thickness of the epitaxial layer. Since the assumption 

of an abrupt junction implies a constant donor concentration, P (x) is 

n 

constant, and thus 
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R 


V 


- '<0 
A 


(. 2 . 12 ) 


Equation 2.12 can be expressed In teras of bias voltage explicitly. 
Using Eqs. 2.3 c'nd 2.4, the following equation i:i obtained 


R 


▼ 




(1 + 


(2.13) 


The total series resistance can now be expressed in the fora of 

R^ = - R2(1 + (2.14) 

where R^ = ‘^^2 “ ^n'^^'"-'0‘ 

For a real varactor, Sq. 2.10 can oe evaluated nwierically if the 

diff'jsion process is known. However, for the purpose of characterization, 
R^ can be siaply expressed as 

Rg = R^ - R2(1 + (2.15) 

and unknown paraneters can be deteralned experiaentaHy . 

2.2.3 Lead Inductan ce. All the metallic portions of the package 
contribute parasitic ind^‘i''ance which appears to be in series with the 
junction capacitance. The most significant contribution undoubtedly 
comes from the lead wires or bonding straps which connect the semicon- 
ductor die to one of the end caps, beca'use of their very small cross- 
sectional dimensions. Typically, the connection consists of a single 
place of gold wire with a diameter of 25 micrometers ('-m), or an o’^tho- 
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gonal padr of 25 (im thick by 75 wide gold straps, vfith the center 
attached to the semiconductor die aind both dens to the and cap (ref or co 
Fig. 2.1 a) . 

The inductance of a round wire with len^h 2 and diameter d is given 


as Cl8]] 


L = 



• - 0.75] 


(2.16) 


where L is in nanohemies (oH), and 2 and d are in centimeters. The low 
frequency inductance of a straight rectangular bar with length 2 , width 
w, and thickness t is L13] 


L = Ei[l n ^ + 0.5 + 0.2235 " ■ (2.17) 

where L is again in nH, and all dimension are in centimeters. The induc- 
tance values at microwave freqaencies c *e affected by skin effect, but 
are lower than those given by the low frequency formula by less than 
6 per cent [l9j. 

The value of lead inductance associated with a particular type of 
package can be obtained from the manrf acturer . However, the lead wires 
or bonding straps are almost universally installed by hand, thus the 
lead inductance ijsually 'varies from unit to -unit. Therefore the data 
supplied by the ainiifacturer must be verified experimentally. 

2.2.4 Package Capacitance . In general, capacitance exists between 
any two separated conductors. Therefore, the package capacitance of a 
varactor diode comes from many sources; ‘oe^ween the two end caps, between 
t.he bonding straps and the bottom end cap, between the bonding straps and 
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the mounting post, and between the top end cap and the mounting post. 

Sven if the fringing field effect and discontinuities are neglected, 
evaluation of the package capacitance is still a formidable, if not impos- 
sible, task However, the value of the package capacitance can be 

^ easily obtained by measuring the capacitance of a dimimy package at low 

frequencies. 

I 

I 2.3 LuaT?ed-Sleaent Squl'/alent Circuits 

As depicted in Fig. 2.3(a), the parasitic capacitance and Inductance 
are in fact distributed elements. Thus, the general approach to the pack- 
age representation would seem to be either a three-port transformation 
matrix L2ll, or a distributed-element circuit, as illustrated In Figs. 

2 . 3 (b) and 2 . 3 (c), respectively. However, the transformation matrix 
suffers from two major drawbacks: one, that the matrix elements may not 
have any physical meaning and two, that each matrix evaluation is only 
valid at one frequency. While the distributed-element circuit does 
relate its elements to the physical parameters of the package str'occure, 
it also makes ihe analysis of ihe circuit much more difficult, and t.hus 
renders Itselt 'undesirable. However, if the oackage dimensions are m’uch , 

i 

smaller than one wavelength at the frequency of interest, these diflclen- 
cles can be readily removed by employing a lumped-element equivalent 
circuit. 

To see how the equivalent circuit in Fig. 2.ii{a) is concei'/ed, a 
closer lock at ine distributed-element circuin in Fig. 2.3(c) is in order. 

f 
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Pacisa^d Varactor Hepcresenta:ion 
{ 3 .) parasltics associated with package 
vb, basic three-port representation 

iistributed-elenent equivalent circuit 







The elements and represent the junction capacitance and the 'bulk 
resistance of the semiconductor oiaterial (i. e., the substrate and the 
undapleted epitaxial layer) . The resistance and inductances with su'bscript 
1 through 4 are those contri'buted from the top end cap, the 'bottom end cap, 
the bonding straps, and the mounting post, in that order. is attributed 
to the capacitance between the two end caps, while C 2 is that 'between the 
bonding straps and the bottom end cap, with being that between tie 
bonding straps and the mounting post. 

In Fig. 2.4(a), it is obvious that C . remains unchanged while H 

J s 

becomes the sum of ?/ and all distributed resistive elements. The capac- 
itance between the straps atnd the post in the immediate vicinity of the 
semiconductor wafer forms the first shunt’ element, '<l'th portion of 

the strap inductance and all the post inductance appearing Immediately 
in series with the semiconductor elements past the Initial shunt capac- 
itance as represented 'by T^ie second shunt element, represents 

t^e other part of the capacitance between the straps and the post, toge- 
ther with portion of the capacitance 'between the two end caps in the air 
volume. the second series element, consists of the other part of 

the strap indiictance and all the end cap Inductance. Finally, repre- 
sents the capacitance between the end caps in the ceramic casing. 

This equivalent circuit does, to some extent, actually relate the 


circuit elements to the physical parameters to the package structure, and 
therefore is relatively invariant over a wide fre"'uency range. Thus, in 
man/ applications it is 'useful over several octaves for freq'uencies as 


-24- 


high that In the lower end of the miUinie ter- wave range. For lower 
freq.uenclea, this eq.ui valent circuit can be reduced to that of Fig. 2.4(b). 
This Is done by observing that if two adjacent circuit elements, a series 
inductance and a shunt capacitance, h3.ving small immittances, but not 
T enough to be neglected outrif^t, they may often be interchanged in 
position. This will allow them to be combined with elements of the same 
typo and thus reduce the number of loops in the circuit by one. Consi- 
dering the L-section containing and ^'^ 2 * ABCD matrix (see Appendix 
a) which represents this section is 


-^p2 


1 

0 


J*)L 


s2 


>^2 


^ S2 


( 2 . 18 ) 


If the positions of and interchanged, then the 

matrix becomes 


.\3CII 


- 


- 



1 >> l ;2 


1 0 


X - 

0 1 




1 


(2.19) 


Equations 2.18 and 2.19 indicate that the positions of and 
may bo interchanged without introducing significant error for frequencies 
satisfying the condition 

( 2 . 20 ) 




For most standard packages, the criterion of Sq. 2.20 can be easily 
met for frequencies in the X-baTiC (3-12.4 GHz) or lower. The eq-ji'/alent 
circul: can bo further simplified to these in Figs. 2.4(c) and i) . 
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Because of their slnpllcltles, these two circuits, la one fora or the 
other, have "been used hy numerous authors. However, they are only valid 
at ultra-high frequencies or lower microwave frequencies, and any attempt 
to use them at high frequencies might cause significant discrepancies 
between theoretical and actual responses. 


2.4 Measurement of Packard Varactor Diodes 


The measurement techniques of microwave varactor diodes can be 

broadly characterised into two catagories: the "transmission loss versus 

frequency" method due to DeLoach C22^, with extensions by Roberts and 

Wilson L23J» and the "impedance versus bias" method due to Moulding [242, 

with extensions by Harrison [252* Diode parameters, namely C . and R , 

J s 

are calculated from the measured data and expressed in terms of Q factor 
[262. from which amplifier performances can be roughly predicted. Para- 
sitic elements are not of great importance, and thus are not accurately 
determined, since they can be easily tuned out in waveguide or coaxial 
line circ'Jits. Unfortunately, for the current study, this is not the 
case because of the lack of tunabdllty in microwave Integrated circuits. 

In the following sections, a measxirement techniqxie is described which 
will enable the determination of the diode equl’/alent circuit more accu- 
rately. The diodes used in this study are Microwave Associates' MA 48505E 
gallium-arsenide (CaAs) varactor diodes, commonly employed in microwave 
parametric amplifiers. Figure 2.5 shows its dimensions. Diode parameters 
at two reverse bias voltages, 0 and 6 volts, were supplied by the maru- 
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at 0 volt. 


2,4,1 Measurement Setup . Since parametric amplifiers are generally 
used for low-level reception with signal level usually below -70 dBm 
(10* the measurement must be carried out with power level compatible 

to the '.ow-level condition. Preliminary investigations indicated that 
the dic-ie parameters are affected by power levels when the Incident 
power ’..s above -15 dBm. As a resiHt, the power level under which the 
diodes are to be measured was then decided to be about -20 dBm. This 
requires a aeasiirement system with very high sensitivity, and conseq-uently 


facturer, and are listed in Table 2.1. 


Table 2.1 Manufacturer Supplied Data for MA 48509E Varactor Diode 


Parameter 

Diode #1 

Diode #2 

Diode #3 

Cjo (PP) 

0.578 

0.552 

0.594 

(I*') 

0.271 

0.252 

0.280 

Lg (nH) 

0.3 

0.3 

0.3 

Gp (pF) 

0,292 

0.292 

0.292 

f *(GHz) 

292 

294 

267 

Vj (volts) 

13 

18 

18 










IM, : 
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the superheterodyne system shown In Fig. 2.6 was chosen. 

A Irief description of the system follows. To eliminate the pulling 
effect, a circulator with one port terminated In a 50-ohm load Is inserted 
between the generator and the rest of the system. Frequency Is 

measured by a precision frequency meter. Power level is continuously 
monitored throu^^ a directional coupler. Bias voltage Is applied via a 
coaxial blas-^ee. The lecision slotted- line is a 50-ohm HP type 816 a 
with APC-7 connectors. Signal picked up by the untuned RF pro'oe is fed 
to the low-noise mixer. The IF output from the mixer is then displayed 
by a precisely calibrated, 30 MHz amplifier (GR type 1236) . 

The diode test mount deserves a more detailed description since it 

/ 

plays the most important role in the whole system, and is not commercially 
available. To be compatible with the slotted- line, the diode test mount 
was designed around an APC-7 air line connector. Every part was machined 
in brass and then gold plated to minimize ohmic xoss. The diode is held 
between the inner conductor and a cylindrical slug with the center of one 
end slightly recessed to ensure that the diode will be properly centered. 

A fine thread screw is used to keep the slug, and thus the diode, firmly 
in position. Figure 2.7 is a photograph of the completely disassembled 


^1 


diode mount. Some critical dimensions are indicated in Fig. 2.3. 

2.4.2 System Calibration. Before the diode measurement can be made, 


two important parameters of the system must be determined: the attenuation 
constant of the slotted- line , and the position of the reference plane 


consistent with tne scale on the slotted- line . Losses on a slotted-llne 
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Photograph of the Disassembled Diode Test Mount 
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Fig. 2.0 Diode Test Mount 
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are usually very small, and can be neglected In most cases. However, in 

this case, the diodes Iiave a very small resistance, and hence the voltage 

standing-wave ratio (VSWH) is very high. a result, neglecting the 

line loss can affect the measured data drastically. Sandler L2?]] has 

shown that the attenuation constant, a, of a slotted- line is related to 

the short-circuited VSVR, S , by 

0 

a = — nepers/ cm (2.21) 

^0 

where i is the distance between the position of the probe and the reference 
plane In centimeters. The attenuation constant measirred by this method 
has a value of O.OOO^o nepers/cm at 5 GHz and O.OOO 56 nepers/ cm at 6 GHz. 

To establish the load reference plane, the center- recessed slug was 
replaced by another slug with identical dimensions, yet an even and 
smooth stirface, which made firm contact with the inner conductor. The 
positions of standing-wave minima were recorded at intervals of approxi- 
mately 1 GHz from 1 to 18 GHz. These positions were numerioally extra- 
polated by Integer numbers of half -wavelength ’until they coincided, and 
this, of course, was the load reference plane. The value cocsistent 
with the slotted- line scale, was determined to be -1.383 co. The scatter 
through the whole frequency range was rather snail and could well be 
attributed to experimeLtal errors, thus the system was believed to be 
substantially free from major defects. 

For the diode measured in this report, t.he diameter of the larger 
end cap Is 3*05 ami which is Identical to the diameter of the inner 
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conductor. Thus, the load reference plane can "be further moved to the 
end cap-dlelectrlc casing Interface. As will be seen later In this chapter, 
this makes the data reduction procedure much simpler. Since the thickness 
of the end cap is O.O58 cm, the new load reference plane was then moved 
to -1.441 cm. 


2.4.3 Measurement Tectinlque . Due to the comhinatlcn of low power 
and high VSWR, the accuracy of measiiremaat by the "double-minimum'* method 
L28] is rather poor because of the difficulty in precisely measuring the 
power level at locations of standing-wave minima. To overcome this pro- 
blem, the "four-point" method [29] was employed. Referring to Fig. 2.9, 

this method requires the measinr^ments of four positions (x^, x^, and 

/ 

x^^) and the difference 'oe tween two power levels (? dS) whic.h was 'usually 
taken to be about 10 dB. Measurements were made at intervals of about 
200 MHz from 4.8 GHz to 6.2 GHz, and again from 10 to 12 GHz. At each 
frequency, measurements were made at six different bias voltages: 0, 0.5, 
1.0, 2.0, 4.0, and 6.0 volts. At each bias voltage, t probe moved from 
one end of the slotted- line to the other covering all standing-wave 
minima. Usually several hundred sets of data, were taken for each diode. 

The calculation of standing-wave ratio, 3, is quite straightforward. 
Let d^ be the distance between x^ and x^, and between x^ and x^, then 


r 

S = ; 1 


9X?(0.23026P; - 1 


n 7T dp ^ TT d_| 

sin (-j^) - exp(0.23C26p)sin^'-Y“^) 
« ' « 


1/2 


( 2 . 22 ) 
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where Is the :«,velen^h aud P is the difference in dB between the two 

s 

power levels where the four positions were measured. However, Eq.. 2.22 
gives VSWR at the position of minim un (Xq). For lossy slotted- line, 
the value at the load reference plane Is [28l303 

= cothCarctanh( l/S ) - ci] (2.23) 

where Z is the distance between and the load reference plane. If 
S » 1, and cci « 1, Sq. 2.23 lecomes 


^ ” 1/S - Oi 


( 2 . 2 ^^) 


and the load impedance at nhe reference plane is 


1 - jS, tane 

7 =7 ^ 

^o 3^ - jtand 


(2.25) 


where is the characteristic impedance of the slotted-line and 9 is 
the electric an^e between Xq and the load reference plane, and is given 
as 


9 



( 2 . 26 ) 


A computer program was written to process the measured data according 
to Sqs. 2.22, 2.23, and 2.25. During the initial run, all data points 
were processed individually, ana an "average impedance" was oalculated 
at each bias. The computer printout was carefully examined and data 
point with -/al'ue greatly Hffarent from the average val'ue was considered 
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to be "bad point" dua to experimental errors, and thus was discarded. 
After the cleanup, the data were reprocessed. Figure 2.10 shows one 
page of the computer printout on the initial run. 


I 



I 


I 

> 

1 


2,5 Detemination of Equivalent Circuit Parameters 

2.5.1 Diode Test Mount Equivalent Circuit . As has been discussed 
in great detail by Getsinger [jl], the impedance calculated from Sq,. 2.25 
is in fact the impedance of both the diode mount and the diode itself. 
Figure 2.11(b) shews the diode mount equivalent circuit proposed by 
Getsinger. Circuit elements will now be briefly described. 

The first element from the left, C^, is the fringing capacitance 
between the end cap and the outer conductor. The series inductance, L^, 
is the coajdal inductance caused by magnetic fields in the volume bounded 
by the outer conductor as indicated in Fig. 2.11(a). The pi-network 
formed by G^^, L^, and equivalent network for the radial 

line l 32 ^ with lergt.h h, extending from the diode to the diameter of 
the inner conductor. 

Based on the dimensions of the diode and nhe diode mount, these 
parameters were calculated to be (see Appendix 3) 


C, = 0.C62 

pf 

L = O.llS 

nH 

= =r2 = 

0.025 pF 

L, = 0.058 

nH 


To test its '/alidity , the impedance calculated from this equi'/alent 







?1^. 2.10 Computer Printout of Measured Impedances 
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circuit together with the diode parameters in Table 2.1 are compared 
with the measured impedances. Figure 2.12 shows excellent agreement 
between the calculated and measured reactances, and thus substantially 
confirms the validity of the equivalent circuit. 

Once the equivalent circuit of the diode mount is detemained, the 
diode impedance, Z^, can be eaisily obtained from the measured impedance, 
Z^. The diode reactances at two bias voltages are also plotted in 
Fig. 2.12. 

2.5.2 Comnuter-Aided Ontimization of Diode Parameters and Squi-'/alent 


Gircirlt Slements . The problem oi obtaining numerical ’/alues for the equi- 
valent circuit elements from the measured impedance data is usually solved 
by the least- square polynomial approximation method [33 J* Conceptually, 
it is very soraightforwsird. A polynomial function which fits the mea- 
sured data is compared with a second polynomial function which represents 
the impedance of the equivalent circ’uit. Numerical values of the circuit 
elements are solved from a set of simultaneous equations formed by 
equating the coefficients of like terms from these two polynomial func- 
tions. In practice, this method is very tedious and involves enormous 
amount of computational efforts. Worse yet, the equivalent circuit thus 
determined Is only valid at one bias voltage. 

What is really desired is to determine not just the equivalent circuit 
at any particular bias voltage. Instead, it is more desirable to deter- 
mine the values of the parasitic elements and the diode parameters. 


namely, n, C.^, , and as in 2qs. 2.9 and 2 . 15 . Kncwir.g these 



Im(zJ 
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Fig. 2.12 Comparison of Calculated and Measured Reactances 
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parameters, the equivalent circuit at any given bias voltage can be con- 
tructed. Furthermore, these parameters also permit the more accurate 
calculation of pump power level and amplifier perfomance. 

To determine these parameters, a computer-aided optimization tech- 
nique was employed. This technique can be simply explained as follows. 

Initial parameter values are obtained by "educated guess” or any other 
jiidiclal means. Following a certain -^^trategy, or optimization method, 
each parameter is then adjusted (within certain constraints) so as to 
minimize the error between the measirred impedance and those calculated. 

3y doing this successively, the error is finally minimized, and thus 
the "optimum” peLzrameter values are obtained. This is obviously not a 
systematic method, and a certain amount of experience and intuition is 
required. In fact, many factors such as initial parameter values, step 
sizes, optimization method, and error difinltlon, can all a^ffect the 
final outcome. 

A computer program employing a modified "direct search" method [34j 

* 

was written for this study . The error fiinctlon is defined as 
EHEOR 

* “2LI»(2d,cal.) - (^.27) 

f 

I 

* Two subroutine sim-nar to "DIHECT" auid "EXPLCR" In Appendix C were 
’used in this program to perform t.he "direct search". This method 
will be covered in more detail in Chapter 7 . 



I ^ 
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where and are weighting coefficients, and the summation is over all 
frequencies of interest. As is evident from Fig. 2.10, the real part of 
the meastired impedances is less repeatable, and thus less reliable, than 
the imaginary part. Therefore a value of 0.2 or lass was assigned to 
the ratio of w^/w^. Measured data of each diode were usually optimized 
more than ten times, each with slightly different initial value and step 
size, for the equivalent circuit shown in Fig. 2.4(b). The resxilts were 
then compared and the one with the smallest error was retained. Table 2.2 
shows the optimum parameter values for the three diodes listed in Table 2.1. 


Table 2.2 Optimu/i Parameter Values for MA 4850^ Varactor Diodes 


Parameter 

Diode #1 

Diode #2 

Diode #3 

n 

2.332 

2.154 

2.21a 

4^ (volts) 

1.214 

1.147 

1.126 

CjO 

0.581 

0.552 

0.596 

( ohms) 

1.04 

1.03 

1.11 

(ohms) 

0.12 

0.14 

0.15 

L (nH) 

s 

0.317 

0.324 

0.325 

V (pf) 

0.251 

0.251 

0.247 

'-p2 

0.043 

0.046 

0.051 








where and w^ are weighting coefficients, and the simmation is over all 
frequencies of interest. As is evident from Fig. 2.10, the real part of 
the measured impedances is leas repeatable, and thus leas reliable, than 
the imaginary part. Therefore a value of 0.2 or less was assigned to 
the ratio of Measured data of each diode were usually optimized 

more than ten times, each with slightly different initial value and step 
size, for the equivalent circuit shown in Fig. 2.4(b). The results were 
then compared and the one with the sm a llest error was retained. Table 2.2 
shows the optimum parameter values for the three diodes listed in Table 2.1. 


Table 2.2 Optiir.um Parameter Values for MA 48509E Varactor Diodes 


Parameter 

Diode #1 

Diode #2 

Diode #3 

n 

2.332 

2.154 

2.218 

4* ( volts) 

1.214 

1.147 

1.126 

=j0 

0.581 

0.552 

0.596 

( ohms) 

1.04 

1.03 

1.11 

R 2 (ohms) 

0.12 

0.14 

0.15 

Lg (nH) 

0.317 

0.324 

0.325 


0.251 

0.251 

0.24? 


"p2 


0.043 


0.C46 


0.051 



CHAPTER III THEORY OF NEGATIVE RESISTANCE PARAMETRIC AMPLIFIERS 


J.l Introduction 


In the analysis of parametric amplifiers and converters, a set of 
poKer relations originally ueveloped hy Manley and Rowe [ 35 ] provedes a 
fundamental basis. It is convenient for the purpose of illustration to 
consider the general situation as represented by the circuit of Fig. 3-1* 
Two voltage generators at frequencies f^ and f^ together with associated 
series resistances and ideal band-pass filters are placed across a loss- 
less nonlinear reactance. Each filter presents a short-circuit to the 
desired frequency, and an open-circuit to all other frequencies. In 
addition to the two voltage generators, an infinite array of ideal band- 
pass filteta a.nd load resistances are also connected to the nonlinear 
reactance. Tliese filters are tuned to the various sum and difference 
frequencies which will arise because of the nonlinear reactance . The 
equations that relates the power flowing into (positive power) and out 
of (negetive power) the nonlinear reactance are shown by Manley and Rowe 
to be 


n=-M 


+ nf^ 


(3.1) 




n=0 m=-“ 


•n'T ■+• TTj; - 

1 i 


(3.2) 
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where P 


is the average power flowing into the nonlinear reaotance at 


,n 

the frequency mf^ + nf^*-- Equations 3*1 and- 3*?' are often referred as the 
Manley-Rowe relations-. These relations are independent of the shape of 
reactance characteristics and of the driving power level. 

To illustrate the applications of the Hanley-Rowe relations, a par- 
ticular case will now he exaained in which power is permitted to flow at 
only three frequencies. If tie nonlinear reactance is excited at f^ and 
it will generate a third frequency fy The circuit is assumed to 
present an open-circuit to all other frequencies. It is further assumed 
that the power from the voltage generator at f^ is much smaller than that 
from the voltaige generator at .hiich is responsible for driving the non 
linear reactance, and which is usijally called the pump source. If f^ 
is the defference frequency, i. e., f^ = f^ - f^, Sqs. 3*1 3.2 become 



(3.3) 



(3.d) 


Since pump power is supplied to the nonlinear reactance and is 
positive, it follows that P^ and P^ are negative. This means that the 
nonlinear reactance is supplying power to the voltage generator at 
rather than absorbing power from it. Since this power is independent 
of triat supplied by the generator itself, it follows that infinite power 
gain is possible at f^. This is the case of the so-called negetive 
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resistance parametric amplifiers, amd f^ aad are the signatl and the 
idler freq.uency, respectively. 

The nonlirear reactance can he either inductive or capacitive. In 
this study, the treatment is limited to the latter case in which the non- 
linear depletior'. capacitance of a 'varactor described in last chapter is 
employed. More precisely, a varactor has not only a nonlinear capacitance, 
but also a nonlinear rosls'tanco. Engelhrecht C36] has studied circuits 
containing both nonlinear capacitance and resistance for frequency con- 
verters. The time-varying capacitive and resistive elements were assumed 
to be 90° out of phase , i . a . , they are pumped in time-quadreture . This 
arrangement allows the circuit to exhibit non-reciprocal features not 
found in single nonlinear element circuits. Other have also studied the 
cases in which -the two nonlinear elements are pumped in arbitrary phase 
relationship [37l3Sj. However, for high quality varactors, the nonlinear 
effect in series resistance is insignificantly small, and may thus be 
neglected. Throughout the remainder of this report, the series resistance 
of a varactor will be -treated as a variable ( i . e . , dependent upon bias 
voltage) but linear element. 

In the subsequent sections, the behavior of a varactor under the 
influence of a large-amplitude pump voltage will be investigated. This 
is followed by -the formulation of the small- signal immiftance matrix of 
a pumped varactor. The analysis of parametric amplifiers employing a 
simplified equivalent circuit, i. e., a nonlinear cai-jacitance and a 
linear resistance, has 'oeen covered in a num'oer of weil-known books l 37_, 
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C 38 X 393 L‘^X^^ 3 « numerous journal articles C^2X^3X N’o attempt to 
review the vast literatiure will be made here . Instead, generalized ex- 
pressions for power gain and noise figure will be derived for amplifiers 
employing a more realistic eq_ui valent circuit • ^ 

3.2 Pumping of Varactor Diodes 

In a parametric amplifier, let a large-amplitude pumping voltage of 
frequency f be applied across the varactor. 3y Fomrier series expansion, 
the nonlinear capacitance can then be replaced oy a time— varying linear 
capacitance , 

G(t)=2_^G^e 5 (3.5) 

k=-°° 

The purpose of this section is to calculate the Fourier coefficients 
G, , and the required pun.ping power. Ttie discussion is limited to the case 
of voltage pumping, in which the voltage across the nonlinear capacitance 
is assumed to be sinusoidal. The calculation methods for the case of 
current pumping, in which the current through the nonlinear capacitance 
is assumed to be sinusoidal, are simi.lar to that of voltage pumping, and 
are given elsewhere l38X^1J. 

The range of the pumping voltage, v (t), across tiie nonlinear capaci- 
tance is limited by the contact potential, |i , in the forward cias region, 
and by the breakdown voltage, V^, it the 7.everse 'cdas region, j-ne pumping 
is sail to ’oe full if v (t) varies from -? to 7_. In this case, the bias 
voltage is chosen to 'oe 


If the pumping power Is limited » then the field of variation of v (t) is 

P 

narrower, and the pumping is said to he partial « In this case, the 
voltage should he chosen such that the foJJLowing condition is met at 
any instant, 

-♦ s V^(t) s Tj (3.7) 


Since the voltage across the nonlinear 
voltage pumping, v^(t) can he expressed as 

V (t) = V -K V cosa t 
P ° ? P 


capacitance is sinusoitial in 


(3.3) 


where, without loss of generality, the time origin is chosen in such a 
way that the pumping phase is zero. ?rom 3q. 2.9, the time-varying 
capacitance may he written as 


C(t) = 




L1 + (V +7 oosu t) / 
0 n ■Q-'' 


(3.9) 


Letting a = ’'^p/(’^o becomes 


C(t) = 


G 

JlSL 


(l + a cos<ii t)^'^^ 


(3.10) 


where ~ ^jc/(^ The Fourier coefficients, G^, are then 


given as 
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Siace C(t) is aa even function, it folloKs that 





2 r 

~Y~ 1 C(t}- Gos(ku t) dt 
} ^ 




Jo 


Equations 3.10 and 3.12 yield 


( 3 . 11 ) 


( 3 . 12 ) 




cos(k(i»^t) 

(l + a cosu) t)^^^ 
P 



! C 



-JIa. 


coskx ^ 

1 2ir 

J 

(1 + a cosx)^^^ 


/ 


Aquation 3.13 can he evaluated hy either expauision in tenas of hyper- 
geometric functions L4<’X45], or numeirical integration. Tables 3.1 to 
3.^ give the results of the computation of the first four Fourier coef- 
ficients for several different values of n in the range from 2 to 3. 

The maximua gain- bandwidth product of a negative resistance parametric 
amplifier as derived by Xuh and Fukada [9] is given by 


b(10 ioga, + 6) £ 13.64 vf^ (C^/Cq) 


(3.1^) 


where b is the fractional bandwidth, and G, is the transducer gain. F; 


rom 
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Table 3.1 Ratio of First Fourier Coefficient to DC Capacitance 

as a Function of a 






a 

h * 2.0 

n = 2.2 

^ 0 
n = 2.4 

n-= 2.6 

n - 2.8 

n = 3.0 

0,00 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.05 

1.0005 

1.0004 

1.0004 

1.0003 

1.0003 

1.0003 

0.10 

1.0019 

• 

1.0017 

1.0015 

1.0013 

1.0012 

1.0011 

0.15 

1.0043 

1.0033 

1.0034 

1.0030 

1.0028 

1.0025 

0.20 

1,0077 

1.0068 

1.0060 

1.0054 

1.0049 

1.0045 

0.25 

1.0121 

1.0107 

1.0095 

1.0086 

1.0078 

1.0072 

0.30 

^ 1.0178 

1.0156 

1.0139 

1.0126 

1.0114 

1.0105 

0.35 

. 1.0246 

1.0217 

1.0193 

1.0174 

1.0158 

1.0145 

0.40 

1.0329 

1.0290 

1.0258 

1.0232 

1.0211 

1.0193 

0.45 

1.0429 

1.0376 

1.0335 

1.0301 

1.0274 

1.0250 

0.50 

1.0546 

1.0479 

1.0426 

1.0383 

1.0348 

I.O3I8 

0.55 

1.0686 

i.060i 

1.0534 

1.0480 

1.0435 

1.0398 

0.60 

1.0854 

1.0747 

1.0662 

1.0594 

1.0539 

1.0492 

0.65 

1.1055 

1.0921 

1.0816 

1.0731 

1.0662 

1.0604 

0.70 

1.1300 

1.1133 

1.1001 

1.0896 

1.0810 

1.0739 

0.75 

I.I605 

1.1395 

1.1231 

1.1099 

1.0992 

1.0903 

o.ao 

1.1998 

1.1730 

1.1522 

1.1357 

1.1222 

1.1111 

0.85 

1.2531 

1.2181 

1.1911 

1.1698 

1.1525 

1.1383 

0.90 

I.33I8 

1.2839 

1.2474 

1.2187 

1.1957 

1.1769 

0.95 

1.4737 

1.4001 

1.3450 

1.3025 

1.2687 

1.2414 
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Table 3«2 Normalized Second Fourier Coefficient as a Function of a 



a 

n = 2.0 

a = 2.2 

n = 2.4 

n = 2.6 

n = 2.8 

n = 3.0 

0.00 

0.0000 

■ 0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.05 

0.0125 

0 . 0 U 4 

0.0104 

0.0096 

0.0089 

0.0083 

0.10 

0.0251 

0.0228 

0.0209 

0.0193 

0.0179 

0.0167 

0.15 

0.0377 

0.0343 

0.0314 

0.0290 

0.0270 

0.0252 

0.20 

0.0506 

0.0460 

0.0421 

0.0389 

0.0361 

0.0337 

0.25 

0.0636 

0.0579 

0.0530 

0.0489 

0.0454 

0.0424 

0.30 

0.0770 

0.0700 

0.0642 

0.0592 

0.0550 

0.0513 

0.35 

0.0907 

0.0825 

0.0756 

0.0698 

0.0648 

0.0605 

0.40 

0.1049 

0.0954 

0.0874 

0.0807 

0.0749 

0.0699 

0.45 

0.1197 

0.1088 

0.0997 

0.0921 

0.0855 

0.0797 

0.50 

0.1352 

0.1229 

0.1127 

0.1040 

0.0965 

0.0901 

0.55 

0.1516 

0.1378 

0.1263 

0.1166 

0.1082 

0.1009 

0.60 

0.1691 

0.1537 

0.1409 

0.1300 

0.1207 

0.1126 

0.65 

0.1880 

0.1709 

0.1566 

0.1445 

0.1341 

0.1251 

0.70 

0.2088 

0.1898 

0.1739 

0 . 16 C 4 

0.1489 

0.1389 

0.75 

0.2321 

0.2110 

0.1933 

0.1783 

0.1654 

0.1543 

0.30 

0.2590 

0.235^ 

0.2156 

0.1988 

0.1844 

0.1720 

0.85 

0.2912 

0.2647 

0.2424 

0.2235 

0.2072 

0.1932 

0,90 

0.3329 

0.3025 

0.2769 

0.2552 

0.2365 

0.2204 

0.95 

0.3953 

0.3591 

0.3286 

0.3026 

0.2302 

0.2608 


Table 3.3 Normalized Third Fourier Coefficient as a Function of a .1 


IV=ol 


a 

n a 2.0 

n a 2.2 

n a 2.4 

n a 2.6 

n a 2.8 

n a 3.0 

0.00 

0.0000 

0.0000 

O.OUOO 

0.0000 

0.0000 

0.0000 

0.05 

0.0002 

0.0002 

0.0002 

0.0002 

0.0002 

0.0001 

0.10 

0.0009 

0.0008 

0.0007 

0.0007 

0.0006 

0.0006 

0.15 

0.0021 

0.0019 

0.0017 

0.0015 

0.0014 

0.0013 

0.20 

0.0038 

0.0034 

0.0030 

0.C027 

0.0025 

0.0023 

0.25 

0.0061 

0.0053 

0.0048 

0.0043 

0.0039 

0.0036 

0.30 

0.0089 

0.0078 

0.0070 

0.0063 

0.0057 

0.0053 

0.35 

0.0123 

0.0109 

0.0097 

0.0087 

0.0080 

0.0073 

0.40 

0.0165 

0.0145 

0.0129 

0.0117 

0.0106 

0.0097 

0.45 

0.0214 

0.0189 

0.0168 

0.0152 

0.0138 

0.0127 

0.50 

0.0273 

0.0240 

0.0214 

0.0193 

0.0176 

0.0161 

0.55 

0.0342 

0.0302 

0.0269 

0.0243 

0.0221 

0.0202 

0.60 

0.0425 

0.0374 

0.0334 

0.0301 

0.0274 

0.0251 

0.65 

0.0524 

0.0462 

0.0412 

0.0372 

0.0338 

0.0310 

0.70 

0.0645 

0.0568 

0.0507 

0.0457 

0.0415 

O.O38I 

0.75 

0.0793 

0.0699 

0.0624 

0.0562 

0.0511 

0.0468 

0.80 

0.0983 

0.0866 

0.0773 

0.0696 

0.0633 

0.0579 

0.85 

0.1235 

0.1088 

0.0970 

0.0874 

0.0794 

0.0727 

0.90 

0.1598 

0.1408 

0.1255 

0.1130 

0.1026 

0.0938 

0.95 

0.2215 

0.1951 

0.1737 

0.1562 

0.1417 

0.:=294 






I 






m 








' 

.00 

0.0000 

0.0000 

0.0000 

f 

5 0 
w • 

.05 

0.0000 

0.0000 

0.0000 

■ 0, 

.10 

0.0000 

0.0000 

0.0000 

0, 

.15 

0.0001 

0.0001 

0.0001 

0, 

.20 

0.0003 

0.0003 

0.0002 

0. 

.25 

0.0006 

0.0006 

0.0005 

0. 

,30 1 

0.0011 

0.0010 


1 0. 

.35 ' 


0 . 0016 

0.0014 

0. 

i4^ 1 



0.0022 






0. 

45 ( 

D.0042 



0. 

50 ( 


0.0053 ' 


0. 

55 ( 

X ^ 


0.0074 ( 

0.CO65 

0. 

GO o.ous 

0.0102 ( 





0.0220 




0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0002 

0.0002 

0.0002 

0.0004 

0.0004 

0.0004 





0.0011 

0.0010 


0.0017 



0.0026 


0.0041 

0.0037 


0.0058 

0.0052 



0.0072 


0.0109 

0.0098 


0.0149 

0.0134 

0.0121 

0.0203 

0.0182 

0.0165 

0.0278 

0.0250 

0.0227 

0.0390 

0.0350 

0.0317 

0.0569 

0.0510 

0.0462 

0.0913 

0.0813 

0.0740 
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Sq. 3*14, it is o^ous that a large ratio is hi^ily desirable. 

G^/Gq increases, as is evident from Table 3«2, with (l) decreasing n, and 
(2) increasing a. The first condition dictates the choice of a varactor 
with junction being as nearly abrupt as possible. The second condition 
Indicates the varactor should be puaped as hard as possible. However, 
this raises an laportant questions How bard can a varactor be punped? 

In other words, what Is the reasonable value for a? Obviously the ac 
voltage swing can not be allowed to carry all the way to the contact 
potential because of the onset of forward conduction. Such a condition 
introduces shot noise, and thus must be avoided. At the present tine, 
there Is no theoretical analysis on marinun allowable a. Sxperinentally, 
this parameter appears to be from 0.85 to 0.95, depending on the type of 
semiconductor materials and reverse breakdown voltages. For design 
purposes, a value of 0.90 will be assumed. With this assumption, It can 
be observed from these tables that 0(t) can be well approximated by only 
the first three terms of the Fourier series. 

The current flowing through the nonlinear capacitance is L4l], 


i(t) = G(t) 


dt 


Jlfl- 


Tfa 


(1 + a cosWpt) 


V (-« )sinu) t 
P P 


( 3 . 15 ) 


The pumping power then can be calculated from 
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2TT 






If the cutoff freq.uency at a specified Idas V Is defined as 

o 


U s 


c H G., 

S jVg 


(3.17) 


and the noiaallzatlon pumping power at a specified hlas V Is defined as 


(V + 

P =— 2 1 

n R_ 


(3.18) 


then Sq. 3«l6 can he written as 


u 




2rr 


2TT '■ 0) ' 

c 


(a slnx)^ ^ 


/O 


(1 + a cosx) 




(3.19) 


2 H 




2 2ir 


(a slnx)' 


(1 + a cosx)^/^ 


( 3 . 20 ) 


The Integral In Bq. 3.20 can he evaliated hy numerical Integration, 
and. resists are given In Table 3 . 5 . 


3,3 Small Signal Representation of a Nonlinear Capacitance 

The relationships between the current flowing through a nonlinear 
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Table 3*3 Nor m al i zed Fvmiping Power ae a Function of a 




a 

n = 2.0 

a » 2.2 

n = 2.4 

n = 2.6 

n a 2.8 
* 

n a 3.0 

0.00 

0.0000 

0.0000 

0.0000 

0.0000 

o'.oooo 

0.0000 

0.05 

0.0013 

0.0013 

0.0013 

0.0013 

0.0013 

0.0013 

0.10 

0.0050 

0.0050 

0.0050 

0.0050 

0.0050 

0.0050 

0.15 

0.0113 

0.0113 

0.0113 

0.0113 

0.0113 

0.0113 

0.20 

0.0202 

0.0202 

0.0202 

0.0201 

0.0201 

0.0201 

0.25 

0.0318 

0.0317 

0.0316 

0.0316 

0.0316 

0.0315 

0.30 ■ 

0.0461 

0.0459 

0.0458 

0.0457 

0.0456 

0.0456 

0.35 

0.0633 

0.0630 

0.0628 

• 0.0626 

0.0625 

0.0623 

O.iMD 

0.0835 

O.O83O 

0.0826 

0.0823 

0.0821 

0.0819 

QA 5 

0.1070 

0.1062 

0.1056 

0.1051 

0.1047 

0.1044 

0.50 

0.1340 

0.1327 

0.1318 

0.1310 

0.1304 

0.1299 

0.55 

0.1648 

0.1629 

0.1615 

0.1603 

0.1593 

0.1585 

0.60 

0.2000 

0.1972 

0.1950 

0.1932 

O.I9I8 

0.1907 

0.65 

0.2401 

0.2359 

0.2327 

0.2302 

0.2281 

0 . 226 ^ 
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capsLcltaiice and the signal voltage across its terminals will now he de- 
veloped. The amplitude of the signal voltage Is assumed to he much 
smaller than that of the pumping voltage. The treatment Is limited to 
the case of negative resistance parametria amplifiers • 1. e., besides 
the pimping voltage, voltages at only two other frequencies are assumed 
to he of slgnlflcaint magnitude. These frequencies are the signal fre- 
quency f , and the Idler frequency f . . The respective voltages and 
s ^ 

currents at these frequencies are 


Vi(t) 


1 ju t 


1 r ^“l^ 


1 r T 

* — L I3 ^ 


2 ■‘■1 ® 


+ Vg e ] 


+ e ^ ] 


+ I e 
s 


♦ ^3 


+ e ^ ] 


(3.21a) 

(3.21h) 

(3.21c) 

(3.21d) 


where the asterisk denotes the complex conjugate. If it is assumed that 
only the first two terms in the Fourier series are significant, then C(t) 
can he expressed as 

Jw t -j“ * 

0(t) = C(, + 0j^( . ^ + e S’ ) 


= o^a + Yj,{ 


Ju t -jw t 

p + e P 


)] 


e 


(3.22) 
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where ~ G^/Gq. The voltage across a noolliiear capacitance and the 
current flowing throu^ it are related by 

i(t) » [G(t)v(t)] (3.23) 

or 

^s ^1 “ It 


Substituting Sq.s. 3«21 and 3*22 into Eq. 3*24, and equating the coeffi- 
cients on the left- and right-hand sides at f^ and f^^, the following 
equation is obtained, 


s 

« 


G« 
s 0 


-JWo 






S 

* 


(3.25) 


In deriving the snail-signal admittance matrix of Eq. 3 . 25 , it is 

implied that all unwanted harmonics are short-circuited. In practice, 

a perfect short-circuit can never be obtained because of the inevitable 

series resistance. A different set of relationships which correspond to 

a condition of open-circuited harmonics cam also be derived by ng 

several additional frequencies into consideration. These additional 

frequencies are * u + u , u, s - w , and u, s 2<u + u . To 
J P s 4 p s 5 p s 

account for the and terms, the second hammonic component must be 
Included in the expression of C(t). Thus 


G(t) = Gq[ 1 + t^Ce P 


P 


^ )] (3.26) 


+ e ^ ) + Yo(® ^ + e 


Hhers Y2 ~ '^® ®®^® procadura as In the derivation of 

Sq.. 3.25* the aaall-slgnal admittance matrix for the case of open- 
circuited harmonics Is obtained as follows; 


» « 

Is 


j“s°0 

>sVo 

^Vl°0 

•>V2°0 

m 


V 

s 


♦ 

II 


-JVi°0 

0 

1 

-JV2°0 

-jVi°o 

0 


* 

''1 


I 3 

s 

jVi“o 


>3Cq 

0 


• 

"3 

(3.27 

♦ 

I 4 



-nvo 

0 

0 

1 

0 




I5 

m 



0 

jVi^o 

0 

0 


''5 

• « 



I and iT can be solved from Eq.. 3*27 by the Perturbation method C37]i 

S A 

and are given as follows; 




J“s=oVl - 


» 4 

V 

s 

“J‘^j_CqY3_(i - Y 2 ) 

m 

m 

• 

♦ 

^1 

m m 


(3.28) 


The small- signal admittance matrix of either Eq,. 3*25 or Eq. 3 >26 
can be used for the analysis of negative resistance parametric amplifiers, 
depending upon whether the unwanted harmonics are more nearly short- 
circuited or open-clrculted. In practice, it is difficult to control 
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thla condition t and thus it can not hs slaul a priori that one is more 
accurate than the other. The difference in amplifier performance calcu- 
lated for both cases will be shown in Chapter V . 

To include the effects of the series resistance and tlie parasitlcs. 
the admittance matrices must be inverted into impedance lua trices. For 
either case the resulting matrix is 


V 


1 

Y 

¥ 

s 

3 





-Y 

-1 




j(^C 


where 


Y = Yn 


G = Gq(1 - 

for the case of short-circuited harmonics t and 


(3.29) 


(3.30) 


Y = Yt( 1 + Y7)(l - Y,) 


G = C„(l - 


(3.31) 


for the case of open-circuited harmonics. 


3.4 Circuit Analysis of Basic Amplifier Configuration 

Figure 3*2 shows a practical configuration of a negative resistance 
parametric amplifier. A circulator is used to separate the input from 
the output since a negative resistance amplifier is in essence a cne-port 
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affiplifier, 1. e., the input and output signals are of the same frequency 
and at the same port. A circulator Is a non-reciprocal device which 
transfers power from port 1 to port Z, from port 2 to port 3* a^d from 
port 3 to port 1, but not vice versa. Of course. It Is possible to 
extract the output power from a load connected In series with the generator 
and the amplifier. This is not desirable, however, since some of the 
output power Is lost In the generator resistance. The gain and noise 
figure of such a configuration are much Inferior to those of an amplifier 
^naif-5 ng use of a clrculator C^l3« 

The boxes labeled f and f . are band-pass filters centered around 

si 

the signal and the Idler frequenc'^ bands, respectively. Pump and bias 

circuits are not shown In Fig. 3*2, since they are assumed to be perfectly 

Isolated , and thus need not be Included in the analysis. Srpresslons 

♦ 

for transducer gatin and noise figure will now be derived. It should be 

pointed cut that these expressions are Intended for computer-aided design 

implementation, thus the step-by-step derivation is presented In such a 

way as to facilitate computer programming. 

3.^.1 Transducer Gain . In practice, the generator resistance E , 

S 

and the load resistance 3^, are always equal to Z^, the characteristic 
impedance of the circulator. Hence, the transducer gain of the amplifier 
is given by 


♦ Transducer gain is defined as the ratio of the power delivered to the 
load, to the available power of the source. 



Pig. 3.3 Equivalent Circuit of a Hegatlve Resistance Paraaetrlo Aapllfler 


-€9- 
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G 





2 


(3.32) 


Nhars Is the input inpedanca looking into the amplifier from port 2 
of the circulator. Ohvioi'aly, for the transducer gain to he computed, 
the input impedance, must first he determined. To accomplish this, 
consider the circuit in Fig.3.3, in which the diode is replaced V its 
equivalent circuit of Fig. 2.4(h). represents the impedance of the 
idler circuit together with the idler hand-pass filter, while the signal 
circuit and the signal hand-pass filter are represented hy an A2CB matrix 
with elements A(ui ), B(iu ), C(ui }, and D(<» ). To facilitate the analysis, 
the nonlinear capacitance is replaced hy'the small-signal impedance matrix 
of Eq. 3>29 with elements denoted as Z^, Z^, ’Z>2it ^d Z 22 * ^ shown in 
Fig. 3 . 4 . 

If Z^ represents the impedance of the circuit to the right of the 
matrix as indicated in Fig. 3*4. then 




= "3 * (3.33) 


rpi 


where the symbol ” | | " is used to designate two impedances connected in 
parallel. The impedance of the nonlinear capacitance loaded hy the idler 
circuit as seen from the signal circuit is 
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^12^21 


Z s Z,, *- ^ 

^ Zoo + Z, 


'22 "1 






( t ) ( I », C ^( Z . - ' I S ") 

“a i 1 jWj^C ' 


(3.3^) 


nhttCtt "tho skstwrlak d/Bao'tes th# coa plo x conjugate. At this point t it is 
perhaps worthwhile to digress from the derivation and *ow why this type 
of aaplifler is called a negative resistance anplifier, and how the am- 
plication is achieved! 

Substituting = Rj_ - f Sq. 3*3^ becomes 


Z„ = —r 


eq 


"sV ^i ■ *^^^1 " 


>3° 


+ {\ - i/\G)‘^] 








(X^ - l/u»j^C)^J 


- l/u^G) 




2i 


(3.35) 


The second term on the right-hand side of 3q. 3-35 o'oviously represents 
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a, negative resistance 






( 3 . 36 ) 


If the gmai circuit is assumed to be lossless and all reactances 
are tuned out, the impedance looking into the s ign al circuit ftom port 2 
is then a pure negative resistance (not necessarily eq.ual to that of 
Sq.. 3.36 since the circuit may contain impedance transformers), 

thus 



(3.37) 


From Eq. 3*32, the transducer gain of the amplifier is 


1 

2 


i 


2c 


(3.3a) 


From Eq. 3.38, it is evident that amplification is achieved. When 

is m a'i ? to equal to Z , the transducer gain becomes infinite. 

c 

Returning to the derivation of the transducer gain, the impedance 
at the varactor terminals seen from the signal circuit is 


Z 

s 


= (- 


J Vp2 






•) 1 1 (H, 



( 3 . 39 ) 


If the elements of the A3CD matrix of the signal circuit axe A(uj^) , 





^ : SI 
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2 (“g)t «(“g)» and input impedance la then 

, _ A(<-^)Z^ H. B(<-J 

^ ^ 'c(a»g)Zg + D('g (3.^0) 

The tranaducer gain can now be calculated from Bq. 3 . 32 . 

The ptreceedln« foimuUtlon is applicable to both nondegenerate and 
degenerate anpliflera. For nondagenerate amplifiers, Z^ la the Impedance 
of an external idler circuit which ia completely unrelated to the 
circuit. However, for degenerate ampliflera, Z^_ la the impedance looking 
into the aignal circuit fiom the terminala of the varactor. Therefore 

D(a)^)Z^ ^ 

i C(u»^^Z^ + (3*^1) 

where A(o)^), 3(u^), G(«j^), and D(uj^) are the ABCL matrix elements of the 
aignal circuit calculated at the idler frequency. 

— S°i” Figure. The fundamental noise in a parametric amplifier 
is due mainly to the thermal noise generated by the sign^T circuit, the 
idler circuit, and the diode aeries resistance. Although numerous authors 
have provided approximate expressions for predicting parametric amplifier 
noise performance [37l38l4ll46l47], for computer-aided design, a more 
general expression for the circuit in Fig. 3.2 must be developed. 

Thermal noise of a resistor results from the random motion of free 
electrons within the resistor. This random motion causes a fluc- 

tuating voltage, or noise voltage, to exist at the terminals of the 
resistor. Nyqulst [43] was able to show, based on fundamental thermo- 
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d yna a l c co n sl d er ations « that the mean aquace noise voltage In the freq,viency 
interval Af generated hy a resistor R at an absolute temperature T is 


e^ a 4kIRdf 

Q 


(3.^2) 


irtiere k is Boltanann's constant (I.38 x 10 “^^ j/^). For analytical 
purpose, a noisy resistor can be replaced by either a noise voltage source 
in series with a noiseless resistor, or a noise current source in 
with a noiseless resistor, as ^own in Fig. 3.5. 

The equivalent circuit in Fig. 3.4 is redrawn in Fig. 3.6 for purpose 
of noise analysis. As Indicated in the figure. Z_ and Z_. are the in- 

1 S 1 X 

pedances looking from into the signal circuit and the idler circuit, 
respectively. 


^Ts * ®Ts '^^i 



Vp2 








where 


I3(Ug)Z^ + B(wJ 
sm - CCa>g;Z~ + A(a.^) 


(3.43) 


(3.^) 


and 


^T1 “ ^1 


= (- 




i pi 


(3.45) 
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As shoHn in Fig. 3.7, the thermal noise generated by resistances, 

Rts' Rri’ ® 3 » caa now be represented by noise voltage sources, e . 

ns' 

*ai' ®nd* ^®P®ctlvely. These noise sources are uncorrelated and, 
in addition, components of e^^ at which excites the signal circuit 

are uncorrelated to the components of e^^ at which excites the idler 
circuit. 

The matrix equation for the circuit in Fig, 3,7 is 


e , + a 
nd ns 


11 3 Ts 


♦ ♦ 

®md * ®ni 2, 


i3M) 


Z,« R + Z„. I* 

22 s "Ti -^ni 


The noise eucrent in the slsnal elrcuit due to the noise sources at 
“s' “nd “ *ns' 


(3^ * R3 * t zjj^) - Z^ 




Thus, the output noise power at oj is 


N ■= 
s 


I Z 
nss c 


^^ad ®ns^ ^22 ^3 '^T^ 

I ^^11 * ^ Ts^(^22 * ^s - ^12 ^21 ^ 


(3.46) 
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Slnce and ara uncorrelated, It follows that 


(•ad • *0.)^ - ““(T»»Ta * 


whaxa 


B * affective amplifier bandwidth (Hz), 
3 sig n al circuit tamperature (*^), 
a dloda temperature (®K). 

Eq.uatlon 3*48 then becomes 


N = 




^ 22 * 2^1 


(2n * * ^Ta)(222 + 8^ ♦ - Z^ 


S i m i la r ly, the output noise power at is 


N, = 




'12 


(^11 * * =>5 - ^12 Zn 


where is the idler circuit temperature. 

The noise figure of the amplifier is given by [37], 


(3.49) 


(3.50) 


F = 




(3.52) 


where ~ 290 ^ is the standard noise tamperature. 

The preceeding formulation is again applicable to both nondegenerate 
and degenerate amplifiers. However, as discussed in Chapter I, degenerate 
amplifiers possess the distinction of having two noise figures, namely, 
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the single>sideband and the douhla-aldeland noise figures. The single- 
sideband noise figure can be calculated directly fron Eq,. 3.52, while 
the double-aideband noise figure is one-half of that calculated from 
Eq. 3«52. 


I 


CHASTER IV PBOEEBnSS OP 2HCBOSTHI? TEANSMISSIOM LUES 





k,l Introduction 

In the process of technological evolution « continuous interaction 
exists hettreen techniques of different fields. Often techniqijss developed 
for one field stimulate and promote the development and progress in others. 
Such is the case of the microwave integrated circuit (MIC). 

Prior to the earlj sixti.es, nearly all microwave systems utilized 
waveguides and coaxial lines. Semiconductor devices such as ''aractors, 
tunnel diodes, and point-contact diodes were used in receivers. Varactor 

harmonic generators were used in a few systems as low power sources. 

/ 

But, for the most part, tubes were still the principal microwave power 
sources. Hence, althou^ the MIC was introduced in the fifties L^9 j, 
progress, both theoretical and experimental, was very slow in the decade 
following its inception. 

In the mid-sixties, with rapid improvements in the low-frequency 
integrated e^'cult (IC) techniques such as epitaxial growth, passivation, 
photolithography, and metal deposition by evaporation or sputtering, 
micrcwave transistors as well as transferred electron devices and ava- 
lanche diodes were developed for solid-state microwave power sources. 
Schcttky-barrler diodes and diodes were also developed for receivers 
and control circuits. With all these solid-state devices available, it 
was only natural that attempts toward circuit miniaturization should be 
made. To achieve miniaturiaation , the MIC was of oomse one of the 
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choices considered. Hot only did advances in low-freq.nency IC techniques 
make these solid-state devices possible , they also made the fabrication 
of olcroHave integrated circuits more precise and more reproducible. 

Out went tha old and crude techniques of glue and razor blades, along 
case the sore sophisticated and more reliable techniques of metal depo- 
sition and jiiotolithogzaphy. "^In ffie years following 19^5* several semi- 
conductor manufacturers, radar system manufacturers, and government lab- 
oratory started programs to develop miniaturized microwave circuits for 
applications in phased-array radar systems. This early work was culmi- 
nated in the .construction of a feasibility radar system by Texas Instru- 
ments in 1968 Tnis radar system uses more than six hundred transai 

receive (T/H) modules fabricated ^ MIC techniques. Each module contains 
a power amplifier, an I? amplifier, a varactor frequency multiplier, a 
balanced Schottky-barrier diode mixer, FIN diode switches, FIN diode 
phase shifters, and associated piiase shift logic circuits. The module 
is about 1.3cm by 6.4 cm, and delivered O.5 at 9 GHz. 

Vita the feasibility demonstrated, nearly all manufacturers of micro 
wave eqiiipment entered the MIG field. Today, MICs are well accepted and 
are employed in essentially all the low- and mediun-power microwave 
applications. 

Microwave integrated circuits can be broadly divided into two cata- 
gories: monolithic and hybrid. A monolithic MIC is one in which all com- 
ponents, active and passive, axe formed simultaneously on a single piece 
of semiconductor material. Dioe to the large varie-oy of active devices, 
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It Is extrefflel 7 difficult, if act impossible, to specify a diffusion pro- 
cess that simultaneously optimizes the performance of every active device^ 
Furthezmore, the interconnections between active devices have relatively 
laxg^ limensiona, especially when distributed elements are used. This 
requires large substrate areas which are very expensive in the case of 
monolithic technology. 

A hybrid MIC is one in which circuit interconnections are formed by 
metal lines on a dielectric substrate and active devices are attached to 
the substrate. Depending upon the techniq.ues by which the metallization 
is formed, the hybrid technology itself can be divided into two classes; 
thin- and thick-film. In thin-film technology, the conductor films are 
deposited in vacuum by evaporation or sputtering, followed by electro- 
plating to increase the film thickness, if necessary. The desired pattern 
is defined by photolithographic techniq.ues. In thick-film technology, 
the metallization is usually formed by silk-screening, in which the 
desired pattern is incorporated into the screen. 

Hybrid MICs can also be classified according to the types of circuit 
elements employed, namely, lumped and distributed* Lumped-element circuits 
for applications at frequencies in the X-band have been reported L5lX32j. 
The limitation on the use of lumped-element clrctdts comes from the ex- 
cessive circuit losses at higher frequencies. 

For the distributed-element circuits, various types of circuit con- 
figurations are currently in use. Figure shows the configurations 
of several commonly used circuits; stripline L33-» microstrip C49_, 


slot 
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liae C34-X553* aai coplanar waveguide [56]. Among these, microstrip la 
certainl/ the most popular configuration, and Is most frequently Iden- 
tified with the term "microwave integrated circuit" . 

In the subaequent sections, electrical characteristics of microstrip 
lines relevant to this study ^fill ^ investigated. These include charac- 
teristic impedance, effective dielectric constant, conductor and dielec- 
tric losses, and frequency dispersion effects. Discontinuities in circuit 
structure, such as open circtiita and T-junction, will be examined in de- 
tail. Finally, analysis and synthesis methods for one circuit component, 
parallel-coupled microstrip band-pass filters, will be developed. 

In HAa-Ung with various techniques in each area, the emphasis is 
placed upon those that are applicable to computer-aided design. In other 
words, closed fora expressions which, in most cases they axe derived 
empirically, are preferred over the more rigorous and more time-consuming 
analytical techniques. However, major works in each area will still be 
referenced, and interested readers are urged to consult the original 
publications. 

Material selections and fabrication techniques, two important aspects 
in MIC technology, will not be covered here; they can be found in books 
dealing with low-frequency ICs or MICs L57J[5SX39 j* 

k,2 Sleetrical Characteristics of Microstrip Lines 

4..2.1 Characteristic Impedance and Effective Dlelecnric Constant . 

The difficulty in microstrip analysis stems from the fact that electro- 
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ffla^etic waves propa^ta alon^ the line in two regions with different 

dielectric constant, i. e., in the substrata with a dielectric constant 

equal to € , and in the air with a dielectric constant equal to 1. In 
r 

such a composite structxsn, the concept of effective dielectric constant 
is quite useful. The effective dielectric constant is a wei^ted mean 
of the dielectric constants of the two regions, and has a value in the 
range from to 1. If the substrate and the air are both replaced by 
a material with a dielectric constant equal to the affective dielectric 
constant, the capacitance between the center strip and the ground plane 
remains unchanged. 

From basic electromagnetic theory, it is obvious that a transverse 
electromagnetic (TSM) wave can not exist in this structure since the wave 
velocity in the substrate is different from that in the air. In fact, 
one can show that not even the pure transverse electric (TE) cr transverse 
magnetic (Th) waves can exist alone in the structure, but that it can only 
support a hybrid mode in which both the longitudinal electric and magnetic 
components are non- zero l60j. The effective dielectric constant therefore 
is dispersive, or frequency-dependent. 

In mlcrostrip analysis, two different approaches are usually taken. 

In the "quasi-ISM” analysis, the structure is assimed to support a TEM 
Wc' /e. The problem then is reduced to one of determining the electrostatic 
potential from a two-dimensional Laplace's equation with proper boimdary 
conditions. Various techniques have been used to solve this problem. 

These include the conformal mapping method by Wheeler finite- 
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element method (relaxation method) by Stinehelfer [,6z2, the moment method 
by Adams [ 63 !] « and the variational method by Yamaslta and Mlttxa [642. 
Dispersion effects are accounted by separate eq.uatlonSt often obtained 
empirically. In the "wave theory" analysis, the dispersion effects are 
calculated directly from a hybrid mode analysis of the structure . Various 
analytical technlq.ues have also been used In this problem [65X662[672. 

Most of the analytical technlq.uas mentioned above require lengthy 
computation that Is generally reserved for a digital computer. No attempt 
to review these techniques will be made here. Instead, analytical and 
synthetic equations by Wheeler [6lj and Schneider [682, with extension 
by Hammerstad [o92< will be introduced. These equations are in closed 
forms, and thus are suitable for computer-aided designs. 

The effective dielectric constant, of a microstrip with width 
W, and substrate thickness H, is given by 

%ff = (4.1) 

for the case of 1, and 

Sff = + 0.04(1 - W/h)^2 

(4.2) 

for the case of W/^ i 1. The characteristic impedance is given as 

z = -==r MaiA + W/^) (4.3) 

^ eff 

for W/H £ 1, a.nd 


(4.4) 


2^ a + 1.393 + 0.667in(wy^ + 1.444)7^ 

9f£ 


for *1. 

For the purposes of syathesls* V/^ can he expressed in terms of 
desired characteristic impedace 2^. For W/fi ^ 2, 


and for 2, 


O' ^ A 

W/H = -|-{3 - 1 - in(2B - l) + — [in(B - l) + 0.39 - 0.6l/€ ]} 

r ' 

(^.6) 


wtisre 2 ^ € - l 

A = -3|o-V2(€^ + 1) ' + -^5 __'o. 23 + O.H/€^j 

r 

and 

3 _ 591.766 
^ Zo 


(^.7) 

(4.3) 


Values of characteristic impedance and affective dielectric constant 
calculated from Eqs.4.1-4.4 are plotted in Fig. '4.2 for several dielectric 
ma-terials, namely, Durold 5380* (^^ = 2.22), fused q.uartz = 3-32), 
and alumina (€ s 9.0), 

2T 


* Trade name for a non-woven glass microfiT»r-reinforced polytetrafluor- 
ethylene (?TFE) structure made hy Rogers Gorp., Chandler, Arizona. 




Effective Dielectric Constant Effective Dielectric Constant 



(d) 

Fig. 4.2 Characteristic Impedance(a) , and Effective Dielec trtc 
Constant(b),(c), and (d) as a Function of W/h. 
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An iogenious empirical equation has been developed hy Schneider 
to describe the dispersion effects. The normalized phase velocity is 
erpcessed as 


where 


V 

P 





(^. 9 ) 


f 


n 



(4.10) 


and and are the microstrip and the £ree-space wavelengths, respec- 
tively. Graphical plots for the normalized phase velocity are shovm in 
Fig. 4 . 3 . It is evident that, for the same W/H ratio, the dispersion 
effect is more profound for a substrate with a higher dielectric constant. 

4.2.2 Conductor and Dielectric Losses . The losses in microstrip 
lines are due to the finite resistivity of the center and ground con- 
ductors as well ais the dissipation in the dielectric substrate. The 
attenuation constant a cam be expressed as [72 ^ 1 


a = 


2P 


(4.11) 


where auid are the powers dissipated in the conductors and the 
dielectric substrate, respectively, and P is the power transmitted along 
the line. Letting a = + a^, then 

■D 


2F 


(4.12) 
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and 


2P 


(4.13) 


If the surface current density distributions on the strip conductor 
and the grovmd plane axe known. Sq,. 4.12 can be written as 


r 


« - ^ 
2Z 

o 


ri(*)| 


/* 


ds + 


‘s2 


2Z 


>(x)| 




dx 


(4.14) 


where and = (irfji^p^)^'^^ denote the surface a *n 

resistivity in ohms per sqaire for the strip conductor and the ground 
plane, respectively, J^(x) and J^(x) the corresponding surface current 
densities, and I the total /surrent per conductor. The quantities 
and P^, p 2 represent the permeability and bull: resistivity of center and 
ground conductors, respectively, and f denotes the operating frequency. 

The first integral is around the periphery of the center conductor, and 
the second integral is over the entire ground plane. Using a technique 
based on the so-called "incremental inductance r'Jle" l 71], Fucel ^ al. 
L?2j were able to derive a set of approximate equations for various line- 
widths. Assuming = R^, the normalised conductor attenuation constant 
is, for V/K S l/2ff , 


2^ H 
c 0 

R 

si 


8.686 

2ir 


•Cl - (W/^)^] [l + H/V‘ + 


H 


V 


+ 4 -)] 


(4.15a) 


for l/2ir 1 W/H 12, 
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^ - (’rf ’M)^] Cl + hA' 


- JL)] 


and for W/k ^ 2, 


(^.15b) 


^ s - 8«-666 ^|- jr ^ tf ■/ H -1 

“si iv'/S. + 0.637fnC2jre(W'/a + 0.9if)]}^ ^ ”wVai + 0.94 


where a = 2 . 71828 • • . . is the case of natural • or Naperian logarithm « and 
t denotes the thickness of the strip conductor. The attenuation constant 
% is in dB/cm. The effective linewidth W is given oy [61], 

W.W+-^( ia-SSL+i) 

for W/H < 1/27T, and 

'''=»♦ 4-( ^=-¥- * 1) (“.i6b) 


iOr i l/2if. i? igura d.d shows plots of the normalized conductor 
attenuation constant versus ¥/1{ for several t/K values, calculated 
using 3q. 1^,15. 

To calculate the die]>.'!Ctric attenuation constant a^, 3q, d.l3 rap 
he shown to equal to 



is 


(^. 17 ) 



-187- 


92. Dell-Imaglne , R, A., "A Parallel Coupled Mlcroetrip Filter Design 
Procedure," G-MTT' 1970 I'nt. Microwave Symp. Digest , pp. 29-32. 

93* Childs, W. H., "Design Techniq.ues for Bandpass Filters Using Sdge- 
Goupled Mlci-ostrip Lines on Fused Silica," G-MTT 1976 Int, Micro- 
wave Symp. Digest , pp. 19 ^ 196 . 

9^. Microwave Engineers' Handhook and Bxiyers* Guide , Horizon House, 
Dedham, Hass., 197^. 

95* Cisco, T. C., Design of Microstrip Components hy Computer , Tech. 
Report NASA-CR-I 982 , March 1972. 

96 . Kaupe, A. F. , Jr., "Algorithm 178i Direct Search," Cbmmun. ACM , 

Vol. 6, pp. 313 - 31^1 June I 963 . 

97. Bell, M., and Pike, M. C., "Remark on 'Algorithm 178: Direct 
Search,'" Commun. ACM, Vol. 9i pp. 684-685, September I 966 . 

98 . Emery, F. E., amd O'Hagen, M., "Optimal Design of Matching Networks 

for Microwave Transistor Amplifiers," Trans, on Microwave 

Theory and Techniciies , Vol. MTT-14, pp. 696 - 698 , December 19^6. 

99 . Bandler, J. N., and Macdonald, P. A., "Optimisation of Microwave 
Networks by Razor Search," LESS Tr^s. on Microwave Theory and 
Techniques , Vol. MTT-17, pp. 552-562, August I 969 . 

100. Denlinger, E. J., et al. , "Microstrip Varactor-Tuned Millimeter- 
Wave IDiPATT Diode Oscillators," LESS Tryis. on Microwave Theory 
and Techniques , Vol. MTT-E3, PP* 953-958, December 1973. 

101. Wu, Y. S., amd Rosenbaum, F. J., '"rfide-Band Operation of Mlcrostrlp 
Ciroulators," IEEE Trans, on Microwave Theory and Techniques , 

Vol. MTT-22, pp. S 49 - 856 , October 1974. 

102. Branner, G. R., and Chan, S.-P., "A New Technique for Synthesis of 

Broad-Band Parametric Amplifiers,’' Trans, on Microwave Theory 

and Techniques , Vol. MTT-21, pp. 437-444, July 1973. 

. Getsinger, W. J., "Coupled Rectangular Bars Between Parallel Plates, 
IRE Tbrans. on Microwave Theory and Techniques , Vol. MTT-10, 
pp. 05 - 72 » January 19^2. 





103 





-91- 


where a is the conductivity of the dielectric substrate. S is the electric 
field intensity at any point inside the substrate, and V is the volta«e 
across the center conductor and the ground plane. The double integral 
is defined over the cross section of the substrate. SLnpson and Tseng 
[ 73 ] have developed an efficient nutterlcal algorithm for evaluating 
Ea. ^<>.17. The calculated results are in exceUent agreement with the 
experimental data of Hyltin M. Table 4.1 gives the normalized dielec- 
tric attenuation constant, a^o. for several dielectric constants. 

Based on the concept of a “filling factor" defined by Poole [?5]. 
Schneider has derived an approximate ea^atlon for dielectric attenuation 

constant L76j, 


0 tan 6 

a, « 27.3 


(4.18) 


where tand = is the loss tangent of the substrate, and is in 

dB per unit length. The filling factor a is given by 


_1 

? - 1 

_ (F + 1] 


(4.19) 


with 


F = (1 + ICH/V)^^^ 


(4.20) 


Numerical results calculated from Eq. 4.18 are also given in Table 4.1. 


Discontinuity Effects in Micro strip Lines 

li.^.l Onen Circuits. In construction of microstrip cir- 
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cuits, an open circuit is usually realised by simply cutting the strip 
off square as illustrated in Fig. ^.5(a). The end region of an open 
circuit formed in this way will store considerably more charge per unit 
length than the remaining portion of the line. - Thus, the end effect can 
be represented by an equivalent capacitance, C , as shown in Fig. ^.5(b). 
This capacitance can be calculated from 

=oc * 

where C(i) is the total capacitance of the section of length Z and width 
W as shown in Fig. b>.5(a), is the capacitance per unit length of a 
uniform line of the same width. In actiial calculation, Z is not infinite, 
instead, it is equal to soma large value beyond which the change in 
C(i) - ZC^ is negligible. 

Alternatively, the end effect can be represented by a length of 

tranaaission line that corresponds to the capacitance C , as shewn in 

00 

Fig. 4-. 5(c). The length hi may be calculated from 

Ai = --g- (i+.22) 

^0 

Many researchers have studied this discontinuity ^77X78X79X80] 
[8lJ, and good agreements are ^s^ially observed between results calculated 
from different techniques, and between calculated and measured results. 

For design purposes, two empirical equations will now be given. 

Silvester and Benedek C78_i have, in addition to the graphical plots, 
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presented their calculated data for opea circuit capacitance by an 
empirical expression 


^oc* 


exp[2O03^Cj^(y( 

k*l 


W xk-1 


H 


)K-X3 


(^.23) 


where the coefficients Cj^ for .six different values of are tabulated 
in Table ^.2. 

Hajaaerstad C69J has found that it Is possible to reduce these equa- 
tions to one valid for all values of by expressing the capacitance 
in terms of an equivalent line extension hi, 


0 419 Vf ^ M/K ^ 0.262 

^ - 0.258 W/H + O.S 13 


IS. 


(h.24) 


Table 4.2 Coefficients c^ for Equation 4.23 
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4.3.2 Microstrip T-jxactlons . The microstrip T-junction, as shown 
in Fig. 4.6(a), is incorporated in a wide variety of microwave circuits. 
The characterization of the tee by a direct jvinction is inadequate since 
it ignores the equivalent reactances associated with the energy stored 
in the nei^hortiood »f the junction. Two different models are usually 
used to represent this discontinuity. In the model shown in Fig. 4.6(b), 
the T-junction is represented by three series inductances and one shunt 
capacitance. These equivalent inductances have been investigated by 
Thomson sued Gopinath CS2J, and the capacitance can be calculated from a 
numerical technique developed by Silvester and Benedek CS3 j. However, 
these computations are extremely time consuming, and thus are not suit- 
able for computer-aided designs. 

The second model as depicted in Fig. 4.6(c) was originally developed 
for stripline T-junction by Franco and 0 liner [843, with extensions by 
Leighton and Milnes C853i Vogel C863, and Hammerstad [69j. The micro- 
strip T-junction is redrawn in Fig, 4.7 to identify the requisite para- 
meters. 

In this model, the most importajit parameter is the stub arm reference 
displacement, i. e., M as indicated in Fig. 4.7. As is evident from 
the figure, Ad may be calculated from 

dd == D^/2 - V^/2 - d^ (4.25) 

where is the li .ewidth of the main line, and d^ are given by 


(4.26) 
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2D, , Z, Z. 

dl * dJo.076 + 0.2(-j^)^ + 0.663 exp( =— i-) - 0.172 

^ X V ^2 ^2 ^2 

(4.27) 

vhttza and Z^ denote the characteristic impedance of the Une and 
the stub, respectively, and effective dielectric constant, 

the wavelength, of the main line. The displacement of the main line 
reference plane is generally very small, and is given as 

= 0.05 n^ D2 Z^/Z2 (4.28) 


where 


Do = 




eff2 ^2 


(4.25) 
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2D- Z^ 
— J — T“ 

^ tsi _2 


^ 2 r- 


z, ^ 


TTT- 

8l 


1 - 


gi 


(4.30) 


The shunt susceptance, B, can be calculated from 


B = 2ir(-^)(2d^/D2 - 


(4.31) 


this modal has been observed to be in fair agreement with experimental 
results L69]. 

In a recent paper by Men2/il and ¥olff [87], the scattering parameters 
of the T-junction are calculated usiig a mode-matching procedure. The 
dy n amic effects are included in t.his model. The discussion on the feasi- 
bility of using this model in computer-aided design will be deferred until 
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Ghapter VI. 

4,4 Mlcrostrlp Para^nal-Goupled Band-Pass Flltera 

Fazallel-coupled mlerostslp Unas have been used extensively in the 

e 

realizations of microwave hand-pass filters and directional couploi's. 

This structure, as shown in Fig. 4.8, can simultaneously -support two 
different modes of propagation, name^.y, even- and odd-modes. In the 
even-mode propagation, waves in the two strips propagate in. the same 
direction; while in the odd-mode propagation, waves propagate in the 
opposite directions. Thus, there are four parameters associated with 
this structure, namely, the even- and odd-mode characteristic impedances, 
and the even- and odd-mode effective dielectric constants. These para- 
meters will be denoted by Z , Z , ^ » and € , respectively, suid 

can be calculated by the analytical techniques developed by Bryant and 
tfeiss C 933 . 

A microstrip parallel-coupled band-pass filter consists of a number 
of .resonators each with a length of approximately half -wavelength . 

Figure 4,9 shows the typical layout of a three-resonator band-pass filter. 
The analytical techniques and design procedure of this type of filters 
will now be given. 

4.4.1 Analysis ’of Parallel-Coupled Band-Pass Filters . To facilitate 
analysis, a single section of parallel-coupled micro strip line Is redrawn 
in Fig. 4.10(a), In which port 1 is connected to the previous section, 


port 3 is connected to the following section, amd ports 2 and 4 are open- 








circuited. Figure 4.10(b) shows the current and voltage associated with 
each port together with the open circuit capacitances of ports 2 and 4. 
The current-voltage relationship of this circuit can be expressed in 
terms of an impedance matrix« 


(4.32) 


where the matrix elements as derived by Zysaan and Johnson CS9j are given 









The even- and odd-mode attenuation constants » <z and a . can he calculated 

Q O 

from Eq.s. ^.15 and 4.17 with replaced hy and Z^^, respectively. 

The even- and odd-mode phase constants« @ and are given hy 

e 0 



2v 

^«o(o) 


(^.35) 


where 1 and 1 are the even- and odd-mode wavelengths, 
fe go 


The current «m.d voltage relations at ports 2 and 4 are 


^2 = ->=cc ''2 


■<f = ->“oc \ 


(4.36a) 

(4.36h) 


Substitution of 3q.. 4.35 into Eq. 4.32 yields (after some lengthy algebraic 
manipulation) , 



where 



tsi 
H • 
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I 

H 

M 

1 
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• 


1 

tsi 

^33 


1 

M 
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(4.37) 


(zf + Z^)(Z^ + 


1 


■itiT 


) - 22^^. 
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For computer-aided analysis, the impedance matrix is then transformed 
to an ABCD matrix (see Appendix A) , 


A 

C 


B 

D 


-I. 


(^.39) 


>ihers 


A = D = Zi^/Zi3 


(4.40a) 


3 = C = 


43 


(4.40h) 


The over-all ABCD matrix for the filter can he readily obtained from 
the individual matrices. Then the input VSWH, Insertion loss, and phase 
shift can be calculated as follows: Assuming both the source impedance 
Z , and load impedance Z^ are real auid positive, the input impedance is 


given by 


“IS CZ^ + D 


7. = 


(it.iti) 


The reflection coefficient is 


r = 




- Z 


4 . * 4 


(lt.42) 


and the input VSVJH is 


S = 


1 + lrl 


( 4 . 43 ) 
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The voltage txansaission coefficient is defined as C90]]* 


T s + jT^ 

= A + B/Z^ + CZg + {kM) 

and the insertion loss and phase shift* in terms of T* are 

L^ = 10fog( |T|^Z/4Zg) (^.45) 

0 = arctan(T^/T^J (4.46) 

where L^, the Insertion loss* is in dB* and 6, the phase shift, is in radian. 

4.4.2 Synthesis of Parallel-Coupled Band-Pass Filters . The design 
procedure was originally developed by Cohn C9U strlpline * and later 

extended by Bell-Imagine C9^j to include microstrip. Cohn has shown that 
each parallel-coupled section with a length of qvarter-wavelength is 
equivalent to an ideal impedance inverter with quarter-wavelength of line 
on either side, ?rom this equivalent circuit, the even- and odd-mode 
characteristic impedances are derived as C 91 j, 



1 + 


o 

K 


* (-K-)' 





2 


(4.4?a) 

(4.47b) 


where is the terminating impedance l^see Fig. 4.9), and X is the imped- 
ance of the ideal impedance inverter, and can be calculated from low- pass 
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prototype aleaeata [9l]I« 

The length of ea^h section ia given hy [92], 


i 



Z 

oe 


+ Z 


00 


"^^effe 


oe 


^V€ 


effo 00 


(4.48) 


This procedure usually results . even for filters with moderate hand- 

width, in very high Z and very low Z for the first and the last sec- 

oe oo 

tions. In terms of physical parameters, this means extremely narrow gaps 
between strips for these sections, and thus requires a hi^ degree of 
accuracy ia photolithographic techniques. However, it has been deter- 
mined that more workable parameters can be obtained by slightly per- 
turbing, using computer-aided optimization techniques, the parameters 
calculated from this design procedure [933 • 

Using this techniqxae, a four-section Butterworth (maxlmally-flat) 
band-pass filter with 2^ bandwidth for operation at 5*5 GHz was designed 
for a Duroid 5SS0 substrate (H = 0.508 mm) . The even- and odd-mode 

characteristic impedances for -each section were first calculated from 

* 

Cohn's equations. A computer program was then used to optimize these 
Impedances. The initial and the optimized values are listed in Table 4.3. 


♦ This program consists of three subroutines similar to "BPFILT", 
"DIHECT", and "EXPLOR in Appendix C. See descriptions of tae.-=je 
subroutines in Chapter V. 
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The width of the first and the last gaps (S^ and In Fig. 4.9) ia in- 
creased from 0.07 mm to 0.11 mm, tl?jis considerably lessening the toler- 
ance requirements imposed on the photolithogra^xic process. The physical 

dimensions are given in Table 4.4. Calculated and actual responses are 

* 

plotted in Fig. 4.11. The passband insertion loss was measured to be 
0.5 dB at 5.5 as compared to the calculated value of 0.38 dB at the 
same frequency. The measured 1-dB bandwidth is 430 MHz (5>30~3.73 GEz) 
as compared to the calculated bandwidth of 460 MHz (5.^~5>73 GHz). 

Table 4.4 also gives the dimensions of a 11 (^z band-pass filter 
with €% bandwidth. Calculated and measured responses of this filter 
were also in excellent agreement. The passband insertion loss was 


approximately 0.7 dB. 
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Table 4,3 Even- and Odd-Mode Characteristic Impedances for 
a Butterworth Band-Pass Filter with Bandwidth 


Optimized Impedances 


Initial Impedances 

Z 

Z 

oe 

oo 

74.01 

38.56 

54,84 

^5.95 


70.15 

55.79 



(AU impedances in ohms) 


Table 4,4 Dimensions of Paxallel-Goupled Band-pass Filters 



W(mm) 

S(mm) 

♦ 

X(mm) 

/ \** 

1.23 

0,11 

9.76 

4.77 

1.51 

0.64 

9.70 

4.74 



* f ^ = 5.5 GHz 


♦* f * U GHz 
0 


















CHAPTER V DESIGN AND REALIZATION OF MIC DEGENERATE AMPLIFIERS 


5.1 Deacriptlog of computer prograa GADDAC 

In recent years, a number of computer-aided design pr ograms have been 
made available to microwave engineers on time-sharing systems [9^1']. ■ Large 
microwave firms have also developed their own. package programs that are 
only available to their design engineers. These p rograms usually require 
large computer memory, and are very costly. Moreover, they. are general 
purpose programs which handle mostly passive networks with some including 
a transistor, and others including a diode as the sole active device. 

The computer program GADDAC (Computer-Aided Design of Degenerate 
Mplifier Circuits) developed in this report is intended primarily for 
the design of degenerate parametric amplifiers, thovigh it is also u"ful 
for certain types of passive networks. It can be readily extended to 
include nondegenerate amplifiers. Furthermore, with a few minor changes, 
the program can bo run on a mini-computer since it requires very little 
memory. The source listing of GADDAC, written in FORTRAN language, is 
given in Appendix G, together with circuit element identification codes 
and Inqut data card requirements. Similar to that of G1scol95j iii struc- 
ture, the program consists of the followings: 

Main pi.'Ogram 
Subroutine HESPON 
Subroitine 3PFTLT 
Subroutine DIRECT 
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MAIM 



Fig. 5.1 Flow Chart of Program CADDAC 















Subroutine EXPLOR 


Subroutine RARIX3H 
Subroutine EVAL 
Subroutine GRAPH 

Subroutine mSPLI (with entry DISP) 

Subroutine OUTPUT 
Subroutine SER5IT 
Functions CSIKKf CCOSH; and CTARH. 

The interrelations between these individual programs are illustrated 
in the flow chart in Pig. 5.1. The arrows are understood to be directed 
from the cal 1 i ng programs toward the called programs. 

Before these programs are explained, a brief review on the basic 
;tructure of a computer-aided design program is perhaps in order. Figure 
5.2 shows a simplified flow chart of such a program. 

The computer reads the circuit topology and initial parameter values, 
and proceeds to calculate the actual response. Next, a comparison is made 
between the actual response and the objective function (desired response) 
from which an error function is generated under certain criteria. At this 
point, the error is tested to sec if it exceeds some prescribed value, if 
it does, the program goes to the optimization subroutine. The optimization 
subroutine generates a set of incremental parameter values, which, when 
added to the previous parameter values, will yield a smaller error. This 
procedure is Iterated until the program is stopped when one of the fol- 
lowing conditions occurs: 



Fig. 5>2 Simplified Flow Chart of a Computer-Aided Design Program 
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(1) that the error is less than a prescrlhed value, 

(2) that a specified number of iterations is exceeded, or 

( 3 ) that an allotted computing time is exceeded. 

Returning to program GADDAC, each subprogram will be now described. 

54 «1 Main program . The main program roads all input data which 
Include objective function, frequencies of Interest, circuit topology, 
circuit parameter values, and values of various parameters associated 
with the optimization procedure. It also decides whether the circuit is 
to be a na lyzed or optimized, whether a graph of gain versus frequency is 
to be plotted, and whether sensitivity analysis is to be performed. 

5 . 1.2 Subroutines HESPON, 3PFLLT. and functions CSm. CCOSH. and 
Q'TANH . This group of subprograms calculate the ^circuit response. The 
three functions calculate valiies of hyperbolic sine, cosine, a-nrf tangent 
functions with complex arguments. 

Subroutine BPFILT calculates the A3CD matrices of the parallel- 
coupled microstrip band-pass fUter at all frequencies of Interest ac- 
cording to the equations derived in Section 4,4.1. To avoid repetitive 
computations, this subroutine is only called by the nia.'i n program when 
the existence of a band-pass filter in the circuit is detected. Once 
the A3CD matrices are calculated, they are stored and made available to 
subroutine RESPON via a COMMON statement. 

Subroutine HESPON calculates transducer power gain according to the 
equations derived in oection 3.^*1* Algorithm for noise figure calcula- 
tions is not included. However, with the equations derived in Section 


3.4.2, the iaplementatlon should be a straightforward matter. 

6.1.3 Subroutines RANDOM. EYAL. DIHECT, andSXPLOR . This group of 
subroutines performs the optimization procedure and thus can be considered 
as the backbone of this program. 

Subroutine RAHIX3H is a random number generator. 

Subroutine EVAL evaluates the error between the objective function 
anti the calculated response. In CADDAC, the error is defined as 

SEECB (3)"'^^/(NFP) (5.1) 


where NFP is the number of frequency points, and S is calculated from 

NFP r- 

® * y . i ^t^^Vcal. 

k=l ^ 

Subroutines IH3ECT and 23CPL0R employ an improved "direct search" 
aexhod to perform circuit optimization. Algorithms L96X97j for the 
original direct search method developed by Hooke and Jeeves L34j are 
available, and modified methods such as “spider search" l98j, "razor 
search" C99J have 2 d.so been published. The algorithm used in CATLAC , 
when compared with others, has been observed to reduce the computing 
time considerably. 

Basically, the direct search consists of two major moves, namely, 
the exploratory move and the pattern move. At the Beg inn i n g, initial 
values are assigned to N circuit parameters which are to be optimized. 
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and error is evaluated at this set of parameter vaXues (called "'basepoint" 
in direct search) . Next an exploratory move is made . The exploratory 
move varies che value of each parameter by some small amount (called 
"step size"), and observes the effect of each of these variations on the 
error function. Those that reduce the error are retained. After the 
ej^oratory move is completed, certain parameters are increeksed, others 
decreased, and still others remain unchan^d. 

A "vector" in the N-dlmensional hyperspace can be defined as the 
difference betneen the new and the initial set of parameter values. A 
move in the direction of this vector is then attempted. This move is 
called a pattern move. If the pattern move is successful in reducing 
the error, a move in the same direction with larger step size is taken, 
and so on, until failure occurs. Upon the failure of a pattern move, 
the last "good" point is established as the new basepoint from which 
another exploratory move is conducted to determine the new direction 
for pattern moves. If this exploratory move also fails, the step size 
is reduced, and the whole process is repeated. The optimization process 
is terminated when the step size is smaller than a prescribed minimum 
step size. The flow charts in Fige. 5.3 and 5*^ outline subroutines 
t}IH£CT andSXPLOR, respectively. 

A number of features that are not shown in the flow charts will now 
be described: 

(l) 3oth the exploratory and the pattern moves are restricted within 
a certain range imposed by parameter constraints. Thi^ Is necessary for, 
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?i«. 5.3 


Flow Chart of Suteoutlne DIrJECT 
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?l 5 . 5.^ Flow Chart cf Subroutine SXPLOR 





otherwise, certain parameters may ac<iuire values that are physically mean- 
losless, such as negative lengths, or values that are physically meaningful 
hut are difficult to he realized in practice. 

(2) An accelerating factor is included in making the pattern moves, 
1. e., each new move*is always larger than the previous move in step size. 

( 3 ) In suhroutlne DIBECT, after an optimization process is terminated 
when the step size smaller than a prescrihed minimum value, it can be 
restarted from a new basepoint generated randomly by subroutine HANIXjM. 

This can be repeated as many times as the user wishes. 

(4) The flow chart in Fig. 5.4 only Indicates the exploratory move 
for one parameter. Actually, each time subroutine SX2L0R is called, all 

N parameters are e^ored at once. This is easily accomplished by a 
DO -loop. 

Subroutine SSN5TT. The effect on transducer power gain by 

the variation of a parameter X is called the gain sensitivity with respect 
to parameter X, , and is given by 



(5.3) 


where x is the nominal value of parameter X. In subroutine SSHSIT, which 

performs sensitivity analysis, the partial differentiation is carried out 
numerically by 




X 


5. 


.l.ers i>x. is taken to 'be 2f6 of the aomioal paraaeter value. The calctilated 
sensitivities are expressed in terms of percentage gain variation with 
respect to one per cent variation in nominal parameter value. 

5 . 1.5 Subrotttines GRAPH. DISPLY, and OUTPUT . This group of subroutines 
does most of the printouts. 

Subroutine GHAIH plots the gain versus freq,uency. response 

Subroutine DISn«f prints out optimized parameter values, number of 
iterations, number of functions evaluated, «uid error. Entry DISF prints 
out new basepoint. 

Subroutine OUTPUT tabulates the gain versus 1'req.uency response. 


;.2 Design Examnles 


Figure 5*5 shows the basic circuit topology chosen for degenerate 
amplifiers to be designed by GADDAC. Quarter-wavelength open-circuited 
transmission lines aure placed in shunt one q.uarter-wavelength behind 
the diode at both the pmp and the signal frequencies. These lines are 
intended to block pump frequency power from reaching the signal output 
port, and the signal frequency power from reaching the pump port. The 
combined insertion losses of these lines together with the respective 
band-pass filters were measured to be nearly 45 d3 at the signal output 


port for pump frequency power, and in excess of 50 at pump port for 



Fig. 5.5 Basic Circuit Topology for MIC Degenerate Amplifiers 
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ffi gna.1 frequency power. They also provide a very high impedance in para- 
llel with the diode at the respective frequencies so that circuit elements 
>lah^nri the diode need not he considered in designing the matching networks. 
This means the «i. gnAl mat(diing network and the pump matching network can 
he designed independently. The pump matching network design can he 
readily accomplished with a Smith Chart, or hy e^qperimental techniques 
[4 i 3. It will not he covered here. 

For gna-i matching network design, program CAI3DAC is used. Although, 
the program was also used successfully in designing amplifiers operated 
at 1.4 CS.Z uaing a Motorola MV 1863D varactor, the discussion la this 
section will he limited to only the 5*5 (Sz amplifiers which employ the 
MA 4850SE varactors listed in Tables 2.1 and 2.2. In the design examples 
to he given below, varactor #2 wan used. Parameters of this varactor at 
1.5 volts bias and 9^ pumping (a = 0.95) axe: 


C . = 0.374 

pF 

R = 0.82 

a 

s 


L = 0.324 
s 

nH 

Ipl = 0.251 

pF 

; - = 0.046 
p2 

pF 

Cq = 0.530 

pF 

\ = 0.367 



The parameter constraints were set as follows; 

Imcedance 15 Cl i Z — 100 Q 

0 

0.1 1 i ^ i 0.5 


Length 
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where 1 Is referred to a freq.uency at the center of the ^gnal frequency 
hand. 

5.2.1 Exaa-ple 1 . In the first example, a simple clrciilt containing 
two series sections of transmission lines and a parallel tuning stub as 
shown in Fig. 5*6 is chosen as the signal matching network. 

The desired response is 15 power gain in the frequency range of 
5.^ Qiz to 5.6 (Slz. Arbitrary initial parameter values of 60 ohms and 
0.3 wavelength were assigned to each element. Table 5>1 shows the 
initial and the optimized parameter values. The calculated frequency 
response is plotted in Fig. 5*7. Power gain of this amplifier fluctuates 
between 1^.3 dB and 15.1 dB. 

/ 

In calculating the frequency response, the condition of short-cir- 
cuited harmonics was assumed for the pumped diode. Using the optimized 
parameter values, an analysis wan mads, assming the condition of open- 
circiiitad harmonics, which yielded a lower gain and a narrower bandwidth. 
The frequency response is also plotted in Fig. 5.7. 


Table 5.1 Initial and Optimized Parameter Values of the Signal 
Matching Network in Fig. 5.^ 


Parameters 

Initial Values 

Optimized Values 

h' h 

60 G, 0.3 1 

25.0 C, 0.2115 A 

^2’ h 

60 a, 0.3 A 

68.1 C, 0.2787 1 

Zy 

60 a, 0.3 X 

25.0 Q, 0.4102 A 








Amplifier Gain 



^ 2.2 Sxaaple 2 . In this example, the signal matching network 
contains two more elements than the previous case, i. e., three series 
transmission line sections and two parallel tuning stubs, as shown in 
Pl«. 5.8. The objective function is 18 dB power gala in the frequency 
range of 5.^ to 5.6 GHz. The initial and the optimized parameter 
values axe given in Table 5.2, and calcuUted frequency response is 

plotted in Fig. 5.9. In the frequency range of interest, gain fluctuates 
between 1? dB and 19,1 dB. 


Table 5.2 Ij^tial and Optimized Parameter Values of the Signal 
Matching Network in Pig. j.S 


Parameters 

Initial Values 

Optimized Values 

Zl. 

55 0.3 1 

52.36 Q, 0.304 X 

^2* h 

55 0, 0.3 X 

36.63 Q, 0.346 X 

Zy 

55 0 , 0.3 X 

52.48 n, 0.441 X 

Z4. \ 

55 0.3 X 

20.11 Q, 0.168 X 

Zy 

55 Q, 0.3 ^ 

53.34 Cl, 0.488 X 


In this example, the sensitivity analysis has been performed. 
Pigure 5.10 shows the gain sensitivity with respect to the matching 
network elements, while Fig. 5 .II shows the gain sensitivities -with 




B’lg. 5.8 Circuit Used In the Second Design Exaaple in Se iMon 5.2 
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XQspect to the diode paiameters. 

In view of these figixres, one obvious conclusion can be Lmedlately’ 
drsim: a degenerate aapUfler is Indeed a very sensitive device. The 
power gain is sensitive to both the circuit and the diode paraaetars. 

This is not lisited to this particular circuit. For the aumerous ampli- 
fier clrcuitsde signed in this study, all with different topologies and 
different objective functions, the gain sensitivities with respect to one 
or more parameters invariably possess some values in excess of 13 . There- 
fore, not only must the diode parameters be precisely measured, but also 
extreme care miist be exercised in the realization of amplifier circuits. 
This point can be easily demonstrated by considering, for example, c\irve 
#5 in Fig. 5.10(b). 

The gain sensitivity with respect to the length of element #5 (see 
Fig. 5.8) is appro:djnately 38 at 5*5 <^Hz where the amplifier has a nominal 
gain of 18.3^ dB. The nominal length of element #3 is 0.488X at 3<5 GHz 
which corresponds to an actual length (on a Duroid substrate with 
1.88 for 2^ = 53*3^ of appro:daately 1.94 cm. Thus, the power gain 
of this amplifier would decrease from 18.34 dB to 11.37 dB if the length 
of element #3 is increased by of 1.94 cm, or 0.194 mm. This should 
serve the purposes of illustrating the importance of precisely determining 
the microstrip discontinuity effects. 

The computing time for the first example in which six parameters 
were optimized, ranged from 20 to 30 seconds with an I® 36 o /?5 computer. 
For the second example, in which ten parameters were optimized, it usually 



Throu^out this stxidyi sevsral dielectric isateriad.s with dielectric 
constants in the low to aediua range were investigated and compared. 
Specifically, they were q.\iartz (6^ = 3,82), polyolefin = 2.32), and 
D\n»id ^ 2.22). Materials with high dielectric constants, such as 
alumiia - 9*0 - 10. 0), were not considered hecause of their more 
profound dispersi'/e characteristics and tighter tolerance requirements. 

Of the three subst ra te materials considered, Duroid has the axhrntage 
over quartz in that the dielectric constant is lower, it is much easier 
to machine and fabricate, and the cost is much lower. However, the lower 
loss of quartz may make it more attractive at higher frequencies. The 
polyolefin material was found to be less stable mechanically than the 
Duroid, and variations in substrate thickness were considerably In excess 
of specified tolerances. 

The amplifier described in Section 5>2.1 was realized on a Duroid 
substrate with a thickness of 0.308 mm. A photograjh of this amplifier, 
housed in a 10 cm by 10 cm aluminum casing, is shown in Fig. 3 >12. Bias 
voltage was applied through a high impedance line (2^ ~ ^0 C) connected 
to the pump circuit. Quarter-wavelength open-circuited transmission lines 
serving as RF chokes, were placed on the bias line one quarter-wavelength 
away from the main line a* both the signal and the pump frequencies. The 
measured insertion losses, irom signal port to bias port at 3o ^Hz, and 
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from pump port to Mas port at 11 GHz, are both In excess of jO dB. 

The pump source employed a 7arian X-13 klystron which can deliver an 
output power of 300 ad^ at 11 GHz. The signal port was connected to a 
Melabs XH-^Zl circulator with an Insertion loss of 0.3 dB from 4.5 GHz 
to 6.5 CRZf and an isolation ranging between 20 dB and 33 dB over the 
same frequency range. Isolation was 24 dB at 5*5 (2Iz. The characteristics 
of the am]iLifier were measured with a Hewlett-Packard Swept Amplitude 
Analyzer (HP 5755A) . 

The amplifier frequency response curve reproduced in Pig. 5*13 wan 
the result of operation with a Mas of 1.5 volts and a pump power level 
of approximately j2 nM at a frequency of 10.995 GHz. Amplifier gain was 
observed to be 1? dB t 1 dB over a frequency range of approximately 70 MHz. 
By increasing the pump power level to 38 » amplifier gain was raised to 

19.5 dst 1 dB over a frequency of 50 before oscillation occurred. 

The frequency agility was demonstrated by successfully operating the 
amplifier over a pump frequency range of 10.7 GHz to 11.14 GHz with Mas 
voltage ranging from O.O 5 to 3*07 volts, and with a minimum gain of 10 dB. 

The amplifier gain characteristics in Pig. 5*13 were recorded at a 
signal power level of -IS dBm (O.OI 6 mH). It was observed ti^t below 
the signal power level of -15 dBm, amplifier gain changed insignif icantly 
with respect to the variations in signal power level, while above -I 5 dBm, 
amplifier gain was decreasing noticeably with increasing signal power, 
and when the signal power was above approximately -9 dBm, amplifier gain 


ceased to exist. 



r 
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The single-sideband noise figure for the amplifier was csilculated, 
for the operating conditions under which the gain curve in Fig. 5*13 was 
recorded, at a diode temperature of 3^5 with a resulting figure of 
3.3 This calculation, however, did not include the transaisslon loss 
in the si gnal filter (about O .5 dB), matching network loss (about 1 dB), 
and circulator loss (about O .3 dB) . Due to the lack of noise source in 
the frequency range of interest, a noise figure measurement was made using 
the simple and crude " signal generator" method l 3?] with a resulting 
figure of 5*7 dB. This method requires that the noise output of the 
amplifier be measured with zero power input and that the power input to 
double the output be measiirad. The output determinations in the presence 
of noise were difficult to accomplished accurately, and hence the accuracy 
of this measuorement was estimated to be 1 1 dB. 

Certain discrepancies betwe«.n the actual and the calculated responses 
were observed. When the amplifier was first constructed using paLrameter 
•/alues from Table 5-l» amplifier gain was observed to be only 7 dB at 5.5 
GHz over a frequency range of about 22 MHz. Additional elements were 
needed to adjust the gain to that shown in Fig. 5 .13. These elements 
(the short stub behind the 5*5 GHz bend- pass filter, and the long stub 
next to the diode as shown in the photograph in Fig. 5.12) were arrived 
at experimentally. This discrepancy may be caused by the non-ideal cha- 
racteristics of the circulator and the microstrip to SMA connecter tran- 
sition. The asymmetrical T-junction near the diode is another possible 



source of trouble . 
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Another discrepancy exists in the bandwidths, namely, the actual 
■bandKldth was only one-third of that predicted. This may be partially 
attributed to the fact that the simple structiire used to separate the 
pump circuit from the signal circuit is only functional over a rather 
narrow freq.uancy range centered at 5*5 outside this frequency range, 
the loading effect of the pump circuit on the signal circuit becomes 
significant, and thus can not be ignored outright. These discrepancies 
will be examined in more detail in Chapter VI together with possible • 
remedies. 


/ 


C3IAPTER VI SUMMSRY, CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


6.1 Siimmary and Conclualoas 

The purpose of this study has been to investigate the featslhility of 
designing aicrowave paraaetric amplifiers by computer-aided optimization 
techniques, with special emphasis placed upon amplifiers in microwave 
integrated circuit form. That such am approach is feasible has been 
demonstrated. The salient features of this study can be sxmimarized as 
follows: 

(1) A precision measurement technique hais been developed for varactor 
characterization. In implementing this technique, a diode test mount 
whicn can accommodate various types of diode packages, has been designed, 
and the test mount equivalent circuit has been accurately determined. 
Measurements have been made of the driving- point impedances of several 

MA ^ 50 ^ '/aractors over a wide frequency range. 

( 2 ) A number of lumped-element equi'/alent circuits for packaged varac- 
tors valid at various frequency ranges have been proposed. Circuit elements 
and diode parameters have been successfully determined from the measured 
impedance data by computer-aided optimization techniques. 

(3) Expressions of power gain and noise figure for parametric ampli- 
fiers employing a more realistic equivalent circuit have been derived. 

These expressions have been presented in such a manner to facilitate 
computer programming. 

(^) Electrical characteristics and discontinuity effects of microstrip 
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transmlsslon lines have been investigated. Analysis and synthesis nethoda 
for parallel-coupled mlcrostrlp band- pass filters have been presented. 

It has been found that the synthesis method usually results in extremely 
high impedances which are difficult to realize in practice. This defi- 
ciency has been removed by employing computer-aided optimization tech- 
niq.ues to perturb the impedances calculated from the synthesis method. 
Band- pass filters with moderate bandwidth have been successfully designed 
and fabricated at 5>5 CIHz and 11 GHz. Actual and calciilatad responses 
have been observed to be in excellent agreement. 

(5) A computer program, GADDAC, has been developed for degenerate 
parametric amplifier designs. It employs an improved "direct search" 
method to perform circuit optimization which has been obser/ed to reduce 
the computing time considerably when compared with other algorithms. 

The program req.clres very little computer memory, and th;is can be modified 
for a mini -computer. 

(6) An MIC degenerate amplifier has been constructed on a Durold 
substrate. The amplifier has a power gain of 17 .5 dS t 1 dB over a fre- 
quency range of 70 MHz. The single- side band noise figure of this ampli- 
fier has been measured to be approximately 5«? dB. 
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a major advanta^ of MIC technology, the integration is desirable, and 
it can be accomplislied by using solid-state devices such as IMFATT diode 
oscillators [lOO], aad mlcrostrip f-Junction circulators [lOl^. Besides 
achieving circuit miniaturization , the integration also serves another 
purpose: it enables one to predict the amplifier response more accurately 
by including the circulator model, either analytical or experimental, in 
the computer program, and thus remove some of the discrepancies between 
the actual and the calculated responses as discussed in Section 

The mlcrostrip T-junction model developed in Section 4.3 is only 
applicable to symmetrical junctions, i. e., the main line Impedances are 
Identical on both sides of the stub. ?or asymmetrical junctions, the 
model of Menzel and Volff [8?J, in which the T-junction is characterized 
in terms of scattering parameters, seems to be adequate. However, for 
computer-aided design purposes, this model suffers from two major draw- 
backs: one, that the calculations of these parameters are relatively time 
consuming, amd two, mhat these parameters are calculated from line widths 
instead of line impedances. This implies that should this model be di- 
rectly implemented in a computer-aided design program, physical parame- 
ters, not electrical parameters from which circuit responses are calculat- 
ed, would be optimized, and thus computing time would be increased dras- 
tically. A possible approach is to first optimize the electrical parame- 
ters by simply neglecting the junction effects. Next, optimized electrical 
parameters are converted to physical parameters, and scattering parameters 
for all junctions in the circuit are calculated at freq'jencies of ir.terest. 


— V 


Optlmlzatioii process is then restarted, but this time only the length of 
each element is to be optimized so that the scattering parameters remain 
valid and need not be calculated over amd over again during the optimi- 
zation process. However, it is not certain that convergence to objec- 
tive function is achievable. 

Although the computer program can perfozm broadband amplifier designs, 
bonsadband amplifier realizations were not attempted in this study. As 
discussed in Section 5*3i the simple structure used to separate the pump 
circuit from the signal circuit is only functional over a rather narrow 
frequency range. ?or broadbeind auapliflers, different structures mi:st 
be used. It is also worthwhile to point out that in the case of broad- 
band amplifier design, the computing time can be reduced considerably 
if the initial parameter values are obtained from one of the conventional 
synthesis methods discussed in Chapter I [1023. 

Finally, on diode characterization , it is highly desirable to have 
the diode measured ^ situ . Intuitively, the lumped-element equivalent 
circuit derived from slotted-line meas’irement with the diode mounted in 
a coaxial test mount should remain valid when the diode is shunt mounted 
in a microstrip if the diameter of t.ne diode package is much smaller 
than the microstrip linewidth. However, if the diameter of the diode 
package is comparable to, or even larger than the microstrip Linewidth, 
the validity of che lumped-element equivalent circuit is doubtful. In 
order to measxnre the diode situ , the launc.hers, 1. e., the ccaxial- 
to-microstrlp transitions, must first be precisely characterized. To 



APPS3JDIX A HEPHESEmTION OF TWO -PORT REWORKS BY ABCD PARAMETERS 


The purpose of this appendix is to give a brief review on ABd) para- 
meters (also referred as general circuit parameters), which are often 
usdd. in two-port network representations. 

In terms of Pig. A.l, the ABCD parameters are defined by the following 
equations 


Vt = AV. + 31^ 

X C 


(A.l) 



(A. 2) 


(A. 3) 



Fig. A.l Definition of Voltages and Currents 
for Two-?crt Networks 





C = slnh(Y2)/Z^ 

y = propagation constant 
= characteristic iapedance 


g. A. 2 A3CD Paraaeters of Seme Comincn Structures 





For certain types of networks, the ABCD parameters are interrelated 
in the following special ways* If the network is reciprocal 

AD - BC = 1 (x.k) 

If the network is symmetrical 

A = D (A.5) 

If the network is lossless, A and D are purely real, and 3 and G are 
purely imaginary. Figure A. 2 gives the A3CD parameters for several 
common structurea. 

These parameters are particularly useful in relating the performance 
of cascaded networks to the performance of each network when operated 
individually. The ABCD parameters of N cascaded networks as shown in 
Fig. A. 3 are given by 



The input Impedance defined in Fig. A.l can be expressed in 
terms of the ABCD parameters and the termination Z^, 


AZ., -r 3 












h. If It is assumed that the fields are not affected by the actual 

configuration of the package, then L is 

c 


where |i is the permeability in air and is equal to 4ir xlO”'^ H/m. 

Substitution of numerical values into 3q. 3. 4 yields 

L = 0.113 nK 
c 


(3 A) 


3 . 3 Evaluation of C , , G , 
== rl-* — 


and L 



The admittance parameters of a radial line pi-network are given by 
Marcuvitz l32j. Special linear combinations of 3essel f’unctlon have been 
defined particularly for this problem, and numerical values are given 
graphically. However, if both h and D, are sma.n in terms of wavelengths, 



ternm 


-15> 


the pi-network elements can be found from a simple coajdal line approxi- 
mation [31] I 





( 3 .. 5 ) 


°r2 ^r * ^ 


•(Df. 

5h" 




(B.6) 


Substitution of numerical values into Eqs. 3.5 and 3.6 yields 

= 0.0577 nK 
= 0.050 pF 

The problem of dividing into G^^ and C^2 ds not a simple matter. 
However, In the present case, this problem may be avoided by obser’/lng 
that at the highest frequency of interest, i. e., 12 GHz, 

u\c ^2)L^ = 0.008 « 1 

Thus the position of and that of either or G^2 interchangeable , 

as discussed previously in Scetion 2.3. It follows that no matter how 

G is divided between G , and C _, the error will always be negligibly 
r ri 

snail. Therefore, it is possible to arbitrarily set 


= C^2 =0-025 pF 


APISNDIX C CADDAC COMPUTES PROGRAM 


G.l Purpose 

This program performs analysis and/or synthesis of lumped-element 
and/or distributed-element degenerate parametric auaplifier circuits. 


C.2 Stnictxire 

This program recognizes ten different structures. These structures 
axe listed below together with their type identification numbers and 
required parameters. All impedances axe in ohms, capacitances in pico- 
farads, inductances in nanohenries, and lengths in fractional wavelengths 
referred to a center frequency in GHz. A negative value for any para- 
meter Indicates that this parameter is to be optimized. However, diode 
parameters, source impedance, and band- pass filter parameters may not be 
optimized. A type number of 0 (zero) indicates the end of network con- 
figiciatlon. 


(l) Series-parallel 3LC 
Type i 1 
Tata 2 R , L , C 



(2) Parallel-parallel HLC 
Type ; 2 
Tata : H, L, C 
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(3) Series-series RLC 
Type I 3 
Data. : R, L, C 

(if) Parallel- series RLC 
Type : 4 
Data : R, L( C 

(3) Series transaiission line 

Type : 3 

Data i Z , I 
o 

(6) Parallel open-circuited 

transmission line 

Type : 6 

Da^a : Z I i 
0 

(?) parallel short-circuited 
transmission Line 
T:rpe I 7 
Data ; Z , i 


(3) P'jmped varactor 


(9) Source inpedance 
Type I 9 


Data 1 R 

8 



(lO) Parallel-coupled 
band- pass filter 
Type 1 10 
Data : See Section C.4 



G.3 Ltaitations 

(1) Number of total, circuit parameters ^ 100 - 

(2) Number of circuit parameters to be optimized S 50 

(3) Number of sections - 50 

(4) Number of frequency points - 101 


C.4 Input Data Caj’ds 


First card 


Variaale 

Columns 

Description 

FCSINTH 

1-10 

Center frequency (GKz) 

FSTART 

11-20 

Starting frequency (GHz) 

FSTCP 

21-30 

Stopping frequency (GHz) 


31-40 

Prequency increment (GHz) 

FFJMP 

41-50 

P'ump frequency (GHz) 


Second card t 


o; 



ZLO 

11-20 

Lower limit of (ohms) 

BHI 

21-30 

Upper limit of length (X) 

HLO 

31-40 

Lower limit of length (X) 

ALPHA 

41-50 

Line loss (dP/X) 

GAINO 

51-60 

Desired power gain (dS) 

Third card : 



DELTA 

1-10 

Initial step size 

DELMIN 

11-20 

■Minimum step size 

DRATIO 

21-30 

Step size reducing factor 

TOL 

31-40 

Minimum error 


Description of circuit topology starts from the fourth card which 
must he the source lapedance card (type 9) • Except for the hand- pass 
filter, each section takes one card and all parameter '/alues are in 
B‘10.0 format. The end of circuit topology description is indicated by 
a type 0 (zero) card. 

For a band- pass filter, a type 10 card is placed in its normal 
position relati'/e to the ot.her sections. Following the type 0 card, 
the parameters of the baind-pass are then gi/on. For a filter with H 
sections, N+1 cards are required. The first card contains: t.he number 
of sections, the conductivity of center conductor ( 10”^aiho/cm) , the 
conducti-rlty of ground plane (lO'^mho/cm) . The second card describes 
the first section, the third card describes the second section, and so 
on. Each card contains: even-mode impedance (ohms), odd- mode Impedance 
(ohms), even-mode effectl-/e dielectric constant, odd-acde effecti'/e 
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dislectrlc coastant, line width (cm), len^h (cm), gap width (mm), open- 

circuit equivalent capacitance (pF). 

Following the filter cards, two more cards are required. The first 

card indicates whether a sensitivity analysis is desired. A letter f on 
• • 

the first column indicates sensitivity analysis is desired, otherwise the 
card is left blank. The second card indicates whether a plot of frequency 
response (usually with smaller frequency increment and wider frequency 
range) is desired. Again, a latter T on t.he first column indicates 
"yes", and blank ot.herwlse. In the first case, FSTART, FSTOP, and FCTIT. 
must also be specified in FIO.O format starting from column 11. 

A set of data cards used for the analysis of the amplifier described 
in Section 5*2.2 are given in-Flg. C.l. 


C .2 


1 1 ,0 

0. m 


Ki.O 


s .s 

5.^ 

b.e. 

C2 

li)0. 

15. 

C . A‘J 

C. lb 

0 . 1 •> 

0.000‘i 

C .2 

0. 1 

. 

bO. 



10. 




■>. 

b2 . 2t6D 

C . 3C 3b 


J . 

3 d . 2 C « 

0 . 3bb 1 


‘1. 


C .bb 1 


t . 

20.11 2b 

C.Hll 


. 

b 3 . 3 3 / ^ 

C .bl77 


. 

bO. 

C . 125 


. 

bO. 

C . 12b 


M . 

O.babd 

C.lttb 

c. a 

0 . 




^.0 

O.bfiHb 



10. I't*) 

0 . H <1 1 

1 . s f b e 

1. 71 Ci 

f<r2 

^6. C« 

1 .bab6 

1 . 792«7 


^ 6 . 1 a 

1 .bfctb 

1 . 7S<.H 

fO . / 4H 

.t,b a 

l . b 6 ? 2 

1 . 7C84 

f s 




1 s 

b.O 

t.c 

C. Cl 


o 

■ri 

TJ 

o 

n 

X 

o 
c : 

r* 






0.3? 

0. ?5l 

0.0^6 



0. 

12 313 

0.‘>db lb 

0. 1 1 19 

O.Ol 755 


111 12 

0.0 to 15 

0.6b37 

0. 02076 

0. 

Ibll? 

0.970 lb 

0.653 

0.02076 

0. 

12127 

0.9HI 3b 

0.1213 

0.0 17 3b 


h'lg. ('.1 Stimple Data Cards i'or Progt'ajn CAULAC 







«-'*** OEGENrSATE PARAMETRIC A.'-PLlFIrRS «*«** 


P5':gRAM =QR C5SIGMNG SIGNAL CI'^CLI'^ CF ?Pi 


S’J = R 

1. 

2 . 


CUT IN = 3 

p rS°C!N 

DIRECT 

EX^LCR 

= ANX'^ 

5VAL 

GPAPH 

cltplt 

D ISPLY 
35NSIT 
ao = lLT 
C ThPEE 


REQUIRED : 


CSINH, CCDSH, A\3 C^i.NH 


''AIN program 

:i''E‘:SirN PAR ( 50 .EL'<ENT( 1C) fG (ICl ) 

CATA RYES/l^'t/ 

CC-'^Zh SIGM 5 0 ,?APva;((5j) ,?Ai?yi.S(50 ) 

PAS Aw( ICC) , IIY?E( 50) ,LIST(5D) 
CCmmcn FPl'^p fFCEN'^R .FSTAR'.FCEL .GAINO 

NPA = AV ,NLI Sf .NLWSEC f NEVA L t N=RE 0 
CC''’' Jn CEL TA tCELMIN ,TCL, dr ATI Cf alpha 

CO ''PL EX A=( ICl) , = F t I CL) ,CF( I Cl) ,DP( 101) 


Q" 




a 



K8P=^3 

READ ( 19 , EGCt =ND=«:':<: » FCENTP , F STA R T ,FS TCP »FDE L » F PU MP 
REAC(rR,5CG) ' fZLC ,0HI .aLCfALPHAtGAINO 
READ ( IP , SCO ) DELTA ,CELVI N tDOATIO, tCL 
5QJ FCP^^AT( 3P13.0 ) 

WRITE! I w, 6 ICC) FSTAR TfFSTCP ,FPLMP 
6100 FCR'^ATt IHU//,2SH PARAMETRIC AMPLIFIER DESIGN,/, 

♦ 15F SIGNAL NETWORK,//, 

♦ 3CH FREOLENCIES CF INTEREST ,/, 

♦ 12F SIGNAL ,F6.3,4H TC ,F8.3,4H GHZ,/, 

♦ 13H PUMPING ,F8.3,4H GHZ,/) 

WRITE! I W, 6200) GA INC ,ALPH A ,ZHI ,ZLC,3HI,3L0 

6200 FCR'MT!17M POWER GAIN ,F3.2,3H D0,/, 

*■ 16H LINE Less ,F7.4'LCH 09/LAM3CA,/, 

^ 26H PARA'4 = tEo CCNSTRAINTS ,/, 

♦ 16H impedance ,/, 

♦ 2CH L?°ER 9CLN0 ,R3.2,5H DHms,/, 

♦ 2CH LCWER 5CLND ,F3.2,5H CH^s,/, 

> 13H length , / , 

*■ 2CH UPPER 3CLN0 ,F3.3,7H LA'-SOA,/, 

♦ 20h lower ECLND' ,F3. 3, 7H LAMBDA,//) 

AL = HA=ALP HA/8 .666 

NFREO=! =STDP-FSTAST )/FDEL+1.6 
10 REAC!IR,5C0) T Y? E , ! ELMENT ! I ) ,I =1 , 6) 

INTYP5=TYP6-*-C.1 
IFIInTYPE) 3C,11C,2C 
20 IF( INTy°E- 10) 5C,fC,2C 

3D Wr ITE! IW.6C0) 

6CD FCR'^ AT( ///, 191- INVALID CARD TYPE , / , I 5 ,3^9 .4 , // ) 

GOTO ^?3 

SO GOTO ! 6C, 6C, 6C , 6C , 7C , 7C ,7C ,90 ,30 , 95) , INTYPE 


I=! EL«EN-! 1 ) ) 62 , 61,62 

£L‘-ENT( 

DC 64 1=1,3 

NPA=iM = ^ip ARAM+ 1 

I^! EL‘^ENT( I ) ) 63,6A,64 
NL IST=NLIST+1 
L I3T(NL IST)=NPARAM 
paRMAX!NLIST)=i:CC.*A3S(EL''SNT(I)) 
PAR’^-ININL 1ST )=D.c:i* A= 3! ELMENT! r ) ) 
Pi = A’-< (‘IP AR AM ) = AE « ( Cl. M =NT ( I ) ) 


ORIGLNAL PAGE IS 
OF POOP 




- 162 - 


70 CO 76 1 = 1,2 

NPiRAy = NP^9AM4. 1 
IF( cLMENtC I ) ) 72,76,76 
72 NL IST=NL IST+l 

L IS7(nL ISDsNPA^AI- 
IF( I .50 . 1 J GCTC 75 
PAPMAXCNL isn = 9Hr 
PAPMlN{NLIST)a3LC 
GOTO 76 

75 parmax(nl rsT)*ZH r 

PARm IN(NL rST)=ZLC 

76 p AR AW ('4PAP AM )=A3 ( C) M=/s^T( J ) J 
^=2 

got:; lc: 

c SCUSCr I'^oecancE CA = C 

I-(\UMS=C+^iPAOA'^+MIST) 61,E2,S1 
"F ITr( I^, 610) 

D AT( //, 26H SCLRCE CARC *'C5T 3E THC crijP7M 

GOTG «=<;a 
NPAS AW = NPAP. AW + I 
PAR AM (NPARAM ) = A3S( 5LWENT( 1) ) 

N=1 

■•'PIT = ( U,62:) 

I FGRMAT(28H INPLT NETWORK CCVF IGUR A TI CN ,// , 

^ 2SM sect T>o5 parameters,/) 

GOTO ICO 

cnoE CAPO 
N = 6 

CC I=1,N 

NPA3 AM=-JPA.RAM.,1 
PAR AM ('IP AR AM ) = A3 S( =L M=MT( I ) ) 

NUm5EC=NlmsEC 
rTYD5(NUMSEC )= INTyoE 

<•' IT=( Iw,630 ) NLMSEC ,INTYPE,(ELmES7(I ) ,t = i 
FC = 'M*(2l4,6P12.<t) - i w 

GC^O 10 

COUPLED LINE BANDPASS FILTER 

NUMSEC = ‘.'LMSEC + 1 
ITYP ECJLMSEC ) = INT'fP = 

K 5P == I 

I”E ( N, 63 5) NLM SEC , IN ■'YPE 
:5 FC -■■!A”(^l4,21)* EANOCA'^ F^LT=R) 


o <-i r> r"» n o o o o 
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C'f poGK 


GCTO 10 

LOO NU.MS5C=NUHS = C + 1 

ITYP = (MLMSEC)*INT>PE 

WRIT=( IW, 620) NUMS5C ,INTyP5 ,(EL^^E^T(I) ,I=1,N) 

GOTO 1C 

no IF(MUMSEC-5C) 12C,13C,12C 
L^O WR ITE( r^,64C) NLMSEC 

640 FCR'4AT(//,29H TQC f*&NY SSCTICNS IN NE TWCRK , 1 5 f // ) 

GOTO <59 3 

120 IF(NP AR VM-ICC ) 15C»15C,14C 
140 WR IT= { IW, 650) NPAR4N 

650 FCR^* AT( //, 20H T2C *<ANY PAR AY = TE= S . 1 5 » // ) 

GC’2 999 

L50 IF(VLIST-53) lt5,1^5,16C 
160 «RrTE( I«,663) NLIST 

660 FCR'-* AT( //, 19H ICC YANY V AR ! AS L = 5 . 1 5 , / /) 

GOTO 99 £ 

PRINT 2UT 9ANCP15' FILTER CCNF IG L'R =T I C\ 

ANO STORE A5CC ''ATPIa CF FILTER AT 
FREQUENCY POINTS CF INTEREST. 

165 IFKEPF ,EQ. C) C-CTO 17 C 
KFILT=0 

CALL 3PFILT(KFILT) 

KF!LT= 1 

RUN A =REQUENCY PESPCNSE 

170 call OLTplT 

Ir(NLIST .ME. C) GOTO 13C 
hR I"E( U. 670) 

670 =CRM AT( ///, L6h END CF ANALYSIS,//) 

GCTO 225 

synthesis required 

150 CC 190 1=1, NLIST 
J = L I3T( I ) 

190 P AR ( I) = PARA.M{ J ) 

CALL DIR. ECT(PAO ,FyER IT) 

IF( IT-TOL ) 2CC ,2CC,21C 
200 w; ITEI IW, 690) 

630 F0:maT( ///, 19H S=A = Cl- SUCCESSFUL) 

GC ' 0 22 0 

21? V,RnE(U, 69 0) 

o90 =:--'*AT(///, :2»- ER=C= A5CVE tclE=ANCE) 

2 20 CDjT'IUE 


r 


QF POOK QUAUTiI 
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WANT A S£NSITI\,IT> AN-ALYSIS? INDICATE Y = S NO. 

R£Ai:(n,513) CHECK 
IF(CHECK-i^YES ) 2AC»23C,2AC 
CALL S£MSIT 

WANT A FREQUENCY RESPCNSS WITH FINER INCBEHENT PLOTTEl 
INDICATE YES OR NC . 

IF YES, SPECIFY- F£TART,FSTCF, AND FOEL. 

REAOdR, 510J CHECK,FSTART,rSTCP,FCEL 
FC3'< ATI AL,9X, 7F1C.C) 

IFICHECK-RYES ) 27C,25C,27C 

N FR EQ = ( -STOP-F START ) / = OELi-l .6 
IF ( K?P= ,EO . : ) GCTQ 255 
CALL cPFILTIKF ILT ) 

FREQ=F5‘^ART 
CC 260 I=1,NF3EQ 
11= I 

CALL c ESPCNI G( I ) ,FREQ ,II ) 

FR EQ= FR cQ ♦"DEL 
NF=NFR=Q 

CALL GRAPH! =STAR T,FCEL ,G ,NF ) 

CCNT INUE 
G 0”0 1 

WP IT = ( IW, 7GC) 

FC = *^ ATI //, 29H OGtTINE TERMINATED BY EPRO = S,//) 

GCTO 1 

STOP 

END 

subroutine = ESFCMGA IN,FR=" ,kSF» 

this subroutine CALCLLATES gain char ACTEP IS’ICS 
CP DEGENERATE PARAMETRIC Amflifie = S 


CIMENSICN 

common 

Cum.mgn 

COMMON 

CCMMCN 

common 

CCM?L EX 

complex 
complex 
CCm?l E < 
COmpl ex 


EETAI2) 

SIGN! 5 0 ,PARM^XI5C) ,PARMIN(5D) 

PAR AMI ICC), I TYPE I 50) ,LISTI5C) 

FPU'M P,F CENT S ,= START, FD5L ,GAING 

NPAF AY ,NLIST,NLmSEC ,NEVAL,NFREQ 

DELTA ,DELM IN ,TCL,DR ATI C, aloha 

Zll ,Z1D .Z21 ,Z22,OI P , ZEC, ZI N ,ZI 

GAMVi ,G S !NM .CGCch ,C "^ANM ,Cmolx ,CCNJG , : J 

AI2),5I2),CI2),D(2),A2I2),C2l2),Z(Z).YlO 

A K 2 ) ,?. H 2) .Cl I 2) ,01 12) , SJI2) ,3J1 12) 

A = ( ICl) ,B = ( LCD ,CF( LCD .CRILOl) 



MMMil 


ORIGINAL PAGE )Z 
OF POOR QUALii Y 


CCMflOM AF,EF,CF ,0 = 


CJ=(0., 1.) 

FIDL£'^»FPU.MP-FR = a 
PI*3.14L5S2654 
0M?GASa2.*P I*FR5Q 
CMEGAI»2.*PI*FICL5R 
BETAt l)=OMcGAS/FCENTS 
BETA( 2 )=CMcCAI/FCFNTR 
SJ 1 1 )=CJ»QM6GAS 
SJ ( 2 ) = cj*':megai 
SJ U 1) = 3J( n^l.E-2 
SJK 2) = SJ(2)*l.E-2 
CO 2 I=lt2 
AU I ) = 1. 
ei( I )=c . 
c 1 (I I = c . 

Cl( I ) = 1 . 

ID0I,NT=1 

nmax=numsec-i 

IF(NMAX-l) 2CC,2CC,5 
CC 130 J = 2,NMAX 
JTY? == ITYO = { J ) 

GOTO ( 10,20,30,^C,5C,eC,7C,13C,13C,aO) , JTyp? 

LU4PEC ... SERIES PARALLEL 



CQ 15 1=1,2 

Z( I )=1 ./PARA>-{ IPCINT + 15"-!. /( SJ( I ) <PAPAM(IPCI NTi-2) ) 
Z( ! )= l./( Z( I )+3J 1( I )*0APA4( IPCINTi-3) ) 

ID0INT= IPO INT43 
got: 150 


L'J'^PED ... parallel parallel 


CC 25 1=1.2 

Yd )=1./PARAM( I?CINT4l»4l./(Sjm *PAPAMnpni.jT + 2) ) 
Y{I) = Y( D+SJK I )»PARAm( iPCIM+3) 

IPGINT=IPGINT43 , 

GCTO 16C 

LJwpEC ... SERIES SERIES 
CO 3 5 1=1,2 

I(I)=RA^A'1( PCIM+D+SJCI )"PA?Ay{ IFCIM4?) 

Z(I) = Z(I)+l./(SJl(I)*PARi''(IPCrMT4.3)) 


(ton ooo OfiO ooo ooo 
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IPlLNT=IPCIM7 + 3 
GOTO 150 

LUMPED ... PARALLEL SERIES 
CO 45 1 = 1,2 

Y( I )»PA.RA,M( IPCir4T+l)+SJ( I )*PARAM( IPCIM>2) 
45 Y( I )= 1 ./( Y(I )♦!./( SJ 1( I )*PAP4f<n PCINT<-3) ) ) 
IPG IMT» IPO IN T + 2 
GOTO 16C 

TPANS.'^ISSION LINE ... SERIES SECTION 
50 CC 55 1=1,2 

ALPHAL=ALPH4«PARAV( IPCINT+2) 

6=tal = 6E‘^a( r )*PlRAf( ICCINT4-2) 

gamm4scmplx( alphal ,betal ) 

A( ! ) = CC0 SH( GAM^<A ) 

C( I )=A( I ) 

B( I ) = CSINr(GAMM4) 

C( I )=9( I )/PA;iAy( IPCINT+1) 

55 0{ I )=B( I )«PA3A.M( IPCINT + L) 

IPO INT= IPG lNT + 2 
GOTO 170 

Tr A*43M ISSION LINE ... SHLNT CFEN 

60 CO 65 1=1,2 

AL3hAL=ALPHA-PARAM( I PCINT-i- 2) 
eE’AL=9ETA(I)»PARAN(IPCINTf2) 

GAMM A=CMPLX[ AL^HAL ,3ETAL ) 

6 5 Y(I )=CTANI-( GAMMA ) /PAR AM IFCIM + 1) 

IP0INT= IPCINT+2 
GOTO 160 

TP A.-jSM ISS ION LINE ... SHLN’ SHC=TE0 
70 CC 75 1=1,2 

ALPHAL= ALPHA«PARAm( IPCINT + 2) 

EETAL = 9ETAn )*OAP AM{ I PC I NT* 2) 
GAMMA=CmPLX(ALPHAL,9ETAL) 

75 Y{I)=l./(C‘^ANI-(3AyMA)*PARAM(IPClNT+l) ) 

IPO INT= I?0 INT + 2 
GC^D 160 

cc'jpleo line bandpass pilter. 

30 :s==<3 = 



ORtGl?;Al. PACr 5^ 
OF POOf{ QL'A'.. lY 


CO 35 1=1,2 
Al I ) = 4r( ISF ) 

Ed ) = 8F( rSF) 

C( I )»CF| ISF) 

C( I ) = 0F{ ISF) 

35 ISF=NFRE0-KSF-^1 
GOTO 170 

4BCD MATRIX OF SERIES ScCTTCN 

150 CO 155 1=1,2 
A( I )=1. 

B{ I) = Z( I) 

C( I ) = 0. 

155 C(n=l. 

GOTO 170 

A3CO MATRIX OF PARALLEL SECTION 

L60 CO 165 1=1,2 
AI I ) = 1. 

Ed ) = 0. 

Cd ) = Y( I ) 

L 65 cm=i. 


MAT= IX '-ULTIPLICATICN 


NETWORKS IN CASCADE 



CO 175 1=1,2 

A2( I ) = A 1( I )*A ( I) *3 1 ( I ) *c ( I ) 
ei(I)=Al(I)=E(I)+ei(I)= 0 {I) 

A l( I )= A2( I ) 

C2( I ) = Ci( I )-l ( I )*01( I )»C( I ) 

CL( I ) = C1( I)»E( I)*D1( I )-D( I) 

CU I ) = C2( I) 

CCiS I i N 'J E 
5 G=? AP AM ( 1 ) 

RL= RG 

Z IC= I 0 1( 2 )«RL+3 II 2 » ) / (C 1( 2) *FL+A 1(2)) 

Z IC= 1 ./Z lO + SJ K 2 )*PAR AM( I PCINT+5 ) 

ZI>1./ZI0 + SJ(2)*FA=’AV( IFCIM-*-4) 

ZIO=l./ZIO + SJl( Zl^PASA'^I I PCI NT>6 ) 

ZIC= l./Z IC<-^AFAM( IPClNT+3) 

Z I3=CCNJG( ZID ) 

AL'X= 1 . -PAR AM I IPCINT + 2 )*PARAV{ IPO I NT<- 2 ) 
Z 11 = 1 ./( 3 J 1 ( 1 ) = ?A = AM(IPCINT+ 1 )=alX) 

Z 12 =PAR:,m( iPOrJT + 2 )/(SJl( 2 ) »pA 5 a m ( i oq + ) » aijx ) 

Z 2 L = -2 I I«PAP AMI !c:I\t+ 2 ) 

Z22=-Z12/PAP AM( l?Crj7 + 2) 
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is 

qOAOITF 


ZEQ=^Ll-^ 12 «i 2 L/(ZIP •^^22 ) 

ZE3=ZH3-*'PAR'M( IPCIMf?) 

ZEQ= l./ZEO + SJ 1 ( I )^P 4 PiV( IPCINTi-o) 

ZEQ=L./ZEQ+SJ( l)*P4Pdf»{ IPCINT+4^) 

ZEQ= l./ZEQ + SJ K I )*OA'»AV( i pcr NT^-5) 

Z£Q»1./Z£Q 

Zr>l=»(Ai(l)<«ZEC+ai( 1) )/(CUL) *ZEQ+C1( 1) ) 
RH0»(ZIM-RG)/(ZIN+RG) 

ARHO»CABS(RHO ) 

GAIN»ARHO*ARHC 
GAIN* IG.«AL0G1C( GAIN ) 

203 RETURN 
END 

SL3^0UTINE goplLTlK ) 

C 

C THIS SU'^PCUTINE CALCULiTES A3CD RA.9A- = Tr = S 2= A 

C CC'JO'.EO LINE 0AJCPASS FILTER AT ALL FREQUENCY P2INTS. 

ZEVEN( 1 C),ZCCC(LC/,WIDTH( 13 ) ,GAOt 2 D) 

K = VEN( L: ) ,KCD0( IC) , LENGTHdOJ ,CCP (10) 
ZL,Z2,Z2,Z4,Z1P ,Z2® ,Z3P,Z4R,X,Y,Z 
A ,3 ,C ,C ,4 1 ,C ItC J.ZTE YF ,C »»PLX 
CS INh,C TANH ,GA«yAG,G AM maE 
SIGN( :C) ,PAR>AX50) ,PARMIN(50) 

°ARAN( ICC) , I TYFE( 50 ) ,LIST( 50 ) 
FRLMP,=CENTS ,FSTART ,FCEL ,GAI NC 
M? 4 ,PAy ,NLI ST ,,NL^ScC ,NEVA L,NrRE 0 
CELTA , 0 ELyiN ,TCL,DFATIC,AL=>HA 
A = { ICl) ,EF( ICl) ,Cr( ICl) , 0 ^( 101 ) 

AF,a = ,CF ,rF 


PI2=6.2S2L65 
CJ= ( 3 1 . ) 

IF( K .N=. 0 ) 30 T2 15 
RcAC( 13,51) R.N,SIGC,SIGG 
M*R.N+0 . 1 

FACTCF*a.CCl*P I2*( SCRT( 1 . /SI GC ) ♦ S CRT ( 1 , /S I GG ) ) 

UFITE( I-^,61) 

ol FCRMAT(///,2ax,2CFFlLTE= CCNF IGUR ET! CN , // , 

+ Alh SECT ZEVEN ZCD0 kEvEN <C0D, 

+ 33I- -»nTH LENGTH GAP ,/ , lOX ,3HCi-'^ , 

♦ 6 X, 3 HCl-‘^, 2 £X, 2 HC^, 7 X, 2 hCV, 1 :x, 2 H'^'^,//) 

C2 13 1=1, N 

3=AC( IR , 51) ZEV = N( I ) ,ZCCD( I ) ,kEVEN( ! ) ,’<:C3 ( I ) 

♦ - I ) , lengtm( I) ,gad( I ) ,c:p( I ) 

51 = A’( 3 = 13 .0 ) 


cr-' e.ns I0N 
real 

CC'^RLEX 
/ C0.NPLEX 

CCNPL EX 
CC^'^CN 
CCH-'^CN 
CCHV0N 
CCNMQN 
CC^NON 
CC»^OLEX 
CCMMGM 
IR=5 
I 6 

VC=29. 57925 


» 
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OF POOK 


URITE(IA,62) I,Z=VEN(I ) tZCDCm t KEVENd ) ,KOOO( ! ) » 
+ '*«IDTH( I ) ♦LE^GTH(I ) ,GAP(I) 

62 FCP..'iAT( I5,2F'?.2.2F':,4 ,p1C.4,F10.3,FI0.A) 

LO C0P( I )»C.C0l*CCP( I ) 

15 FRSQ*F3TART 

CO 30 I=L,NFREC 
AF( I ) = L* 
eF(I)=0. 

CF( r )=o. 

CF( I )= 1. 

0MEGA=PI2«FREC 
RS=FACT0F*SQRT( SPEC ) 

CO 20 J=LfM 

/iLOHAc= 3.5«RS/( w ICTh{ J ) *ZEV5N( J) ) * LENGTH (J) 
iL PH AC= ALP HA E*z EVEN ( J )/ZCDO( J) 

AIJX= 5CRT(KEVEN( J ) ) 

BET AE=C’-'=GA»AU'*LENGTH( J ) / VC 
BETA'1 = 3ETAE»SCP^(KCCC< J» ) /ALX 
GA.HM A S = C?-'P LX I ALP HA E, beta E ) 

GAM'^AC = CyPLX{ ALPHAC.BE'^AC) 

Zl?= ZEVEN( J ) /CTANF(GAV“AE ) 

Z2°= ZQ0D(J ) /CTANt-( GA^MAC ) 

Z3P=ZEVEN(J )/CSINI-(GAy»'AE) 

ZAP= zoocij )/csiNH(GAyyAC) 

Zl=0.5«( Z 1P+Z2P ) 

Z2=0.5*( Z IP-Z 2? ) 

Z3=0.5*( Z3P-Z4P ) 

ZA=0. 5-( Z3°-*-ZAP ) 

ZTEMP^Z Ul. / ( CJ-C*'EGA»CCP( J) ) 

Z=Z3»Z3/ZTEyo-ZT=yP 

X=Z2-Z3-ZA/ZTE'-P 

Y=ZA-Z2*Z3/ZTEyP 

Z1==ZU(Z2«X + ZA«Y)/Z 

Z2P= Z3-*-( Z2*V+ZH*X ) /Z 

A=Z1P/Z2P 

C=A 

B=( Z IP* Z1P-Z2F*ZZP)/Z2P 
C=l./Z2° 

Al- AF( I )*A+?F{ ! )»C 
0P( n = A = ( I ) *3 ♦PF ( I )*D 
AFC I )» a: 

C l=Cr ( I )«A+QF ( I ) *C 
CF( I)=C = ( n*3 ♦OF ( I )*C 
C = l I )=CL 
2 3 CC'i'^IN'jE 

F= E0= FR =C^=0 = L 
CCN'I'' JE 
FE’JCN 


33 


51 

CCMPLEX rUNCTICN CSINl-{Z) 

c 

C hYPERBailC SINS rlNCTICN WITH COMPLEX APGUM5NT. 

COMPLEX ZfCEXP 

CSINH»(C£XP(Z )-CEXP{-ZJ ) /2* 

RETURN 
END 

COMPLEX FUNCTION CCCSH(Z) 

t-YPERcOLIC COSINS FLNCTICN WITH COMPLEX ARGUMENT. 

CCMOLEX Z,C5>o 

CCOSH*(CEXP(Z XEXF(-Z! ) / 2. 

R ETURN 
ENC 

COMPLEX FLNCTICN CTANl-(Z) 

t-YOERSOLlC TANGENT =LNCTICN WITH COMPLEX ARGU»«EN'. 

COMPLEX ZtCMPLX 

X= 3 EAL(Z ) 

Y=AIMAG(Z ) 

I»=((Y.GT. A. 7123) .AND, (Y ,LT. A.712A)) Y-A.7123 
IF(<Y .GT. 1.57C7) .AND. (Y .LT. 1.5709)) Y=1.5707 
CTANH=CMPLX( TAN*^! X ) , TAM Y) ) /CMPLXI 1. , TANH(X) *TAN(Y ) ) 
R ETURfJ 
END 

SU 5 =DUTINE DIP EC T( PAS P ,FMp) 


Co: I'^IZATIDN SLBPQLTINE -- DIRECT SEASCh 


this SJSKC'JTINE P=R=CP'-S CIRCLIT CPTI ''I ZATIOM 
SY JSING DIRECT SEAsch METHOD. 1= THE SEAFCm 
from the original BiSEPDINT REACHES AN 
unsatisfactory result, THIS subroutine ^ill 

GENERATE A NEW 3ASE®CINT RANDOMLY AND RESTART 
the search routine. this will be REPEATED 
MPANOU times. 


01 ‘PENSION PAPP(S:),PAPN{5C)fOARO(50) 
C0‘^‘’0N SIGNI fC) ,PAPWAX{5 0) .PAR'^IMSD) 

COM'^DN para*-! :CC), ITyR=(50) ,LIST(50) 

CO'-'-DN = = U'"P, = CENT= ,FSTA3T,FDcL ,GAINC 
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OF POO'. " 


NP AR A V » N L I S T tNLMScC »%cVAL»NFR5Q 
COMMON CELT A ,OSL'‘IN fTCLfDR ATIC» alpha 

CCMPLEX Ar( lCL)t8F(LCl)»Cr(lCl) *CF(101) 

COMMON AF,RF,CF,0F 

l'^»6 

INITIATION 

MRAN0U»4 
WOITcI !«., 1C) 

LJ FO = MAT( IHl,////, 

♦ 37H OPTIMIZATICN FPCM ORIGINAL RASEPO INT ,// ) 

I X= 77777 

CALL «f anOCMI IX, IV, PNC) 

IX= lY 
<= AN0Li=C 
CEL INT*nELTA 
N£'/AL*0 

2C 20 I=L,NLIST 
23 P AR 3 ( I ) -3AR0( I) 

call EVAL ( FMB ,PAR3 ) 

\EVAL=NEVAL + 1 
NIT= 1 

^0 CC 50 I-'lfNLIST 
SIGN! ' )= 1. 

5 0 ? ARN ( T )-PARP( I ) 

C^LL &V-*L ( Fmn ,PARN ) 

NrV A. ^nEVAL ■*•1 
I=(cmm-T 0L) FCC,5CC,eC 
30 FM?=FMN 

call explcr ( Fm\, parn ) 

IP(=yp_PMN) 4 £C,a;C,IEC 
153 I=(=mn-TCL) 35C,=fC,ieC 
160 FACTCP=1.2 
173 F‘'3==MN 

C 

C PATTERN MOVE 

C 

CC 300 I=l,NLIST 

7CU0*P133( I , 

PAHP ( I )»? ARN( 1 ) 

CELP A = *PARN( I )-TE*<P 

PARM I )*PAS‘i( I )>C£LPAO •FACTOR 

C 

C CHECK =0R RAP AM=7E= ' £ CONSTRAINTS 

r 

I = (RAR‘i( I) -RAO MAXI) ) Z!C,23C,L33 
13: PAR- ( I )=RAR ^AX( I ) 
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GC": 3CC 

220 IF(PAF'H r )-P4PMl\( I ) ) 25C,2CC,302 
250 PARN ( I )»P ARM IN{ I ) 

3C2 CCMTI‘<JUE 

CALL £VAUFMN,PARN) 

N£VAL*N» = VAL +1 
IF(FMN-TOL) 85C»65C,3C5 
3C5 IF(FMN-FMP) 31C,22C,22C 
312 FACTCP«FACTCR tC.2 

IF(FACT2R .GT. 3.) factor-3. 

GOTO I7C 

ic THE pattern *-aVE =AILSf 3ACK UF CN= ST = ? 

AN2 Mi<= -NCTFSR 5>PLC?ATCPY vcvS 

3 22 r'-''l=-MO 

CC •♦00 I-1,NLIST 
A02 PA;n ( I )-PARP ( I ) 

CALL EXPLCR ( PAF'' ) 

IF(=M,\-rnL) a5C,E5C,^2C 
‘♦22 IF(=mp-cm,\) A5C,AfC,UC 

R52UC5 STEP SIZE 

A52 CELTA=0 = LTA-CRATIC 

IS STE= SIZE i’-ALL ENClGHT 

IF( :ELTA-C£LM IN ) 6CCt6CC,5CC 
5 22 CGNTNJS 
N I'-N !’♦ I 
CC'2 9C 

3=I‘.T 21*^ IN’c='^EG lA TE results 

6C2 CC tOA I-l.NLIST 
L = L ISTfl ) 

oCA PA- A>^(L )*PAR9( I ) 

CALL OISPLYf P-o.MT) 

call output 

SEA3CH =RQM PC = vTCLS SET CF EaScPCIMTS haS FAILED 
generate a new set rc 3ASEFCINTS RANOC«LY AND 

re-s’art The seapck r:l”^ine 

TCJ :-^a_rv.p I 6CS,tC=,eC5 

s25 =--*3==-o 

:r 6 26 : = :,nlist 
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606 3i=5l ! >=pab3( I) 

6 09 K.BJi‘jOU=K=ANOU t-I 

I = (<o uNCU— ^‘’■iNCU ) 6lCt6lCf9CC 
610 w= IT = ( I,*, 920) <9AN0L 
920 FC^‘«’A’'( 1H» 

♦ ////»32’r GENERATE A NEv» SET CF BA 3E PCI N'^S ? / t 

^ IgH *»* ATTE'^FT N0.tI2,4H »*»,//) 

CO 6E3 I=1,MLIST 
CALL RANOCMI I>, IT,FNC ) 

IX=IY 

P APF ( I )=0.5*( °AR'<A)a I )+PAPMlN( I )) ♦PAPF( I ) «( 9NC-0.3 ) 
I = ( PARP { I J-PARMa 3(( jj) fc2C,6l5»615 
6L5 PA=P ( T ) = PARMA>( I ) 

GC’O 6 

620 I = ( P ARP ( I )-P AP'^ IN ( I ) ) f25fc25f650 
625 PAR? (U-P ARM INI I ) 

6 50 CCN'FVJE 

rEL"A=2ELr 
CO 660 1=1. aL IS 
L=L 1ST( I ) 

c6 0 =AR A.M(L )=PA=o ( I ) 
k,RI'E( I-.940) 

940 rCPM AT ( //, 17H 'J?V» SA SO PC I NT,/ / , 

23V* SECT I'tPE ■ PARA '‘ETEPSf// ) 

CALL DI3? 

GCTO 40 

550 CO 360 I = l,NLIST 
56 0 PA=? ( I ) = ?AR*1( I ) 

f'AO-s CVI^^ 

CC 370 I=l,'^LISr 
L=LI3^( I ) 

37} P AR AM ( L ) = P A= P ( I ) 

C-C'O 990 

9 CO CC 910 I=l,'iLIST 
I )=PA= 3 ( I ) 

L=L I3T( I ) 

9 iO PAP A^{ L)=PAR p ( I ) 

>PI~E( 1^,920) 

920 FCR'^ AT ( It-I, // , I7h *•» BEST 3ET *••,//) 

fuo- zna 

99 0 CCNTILuC 

CALL OI5PLY( =“P ,M^ ) 

CALL CL’TPlT 
= E" J=‘< 

*** « “V 

SLS-I'JTINE EXPLCF t =•'- 


o o <*> o 
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This SU3SQUTIN5 NAKSS EXPLC5ATCPY f*CV5S 
TC DECIDE ’HE CIRECTICN =CR PAr^EOs mcvES. 


C 


ID 

15 
2 D 
2 5 


3 D 


CI.'^ENS I3N 

CO'^’^CN 

CCHHQN 

CC'^'^ON 

CC' 4 Ma >4 

CC'^HGN 

complex 

CC'^'^ON 


PARM 5C ) 

SIGN! 50 t?AR Y 4 X( 5 C) f PAR'‘IN( 5 «D) 
PAR AM ( ICC 1 , 1 TYPE ! 50 ) ,LI ST( 50 ) 
FPUMP ,PCENTR ,? START ,FOEL »GA! NO 
NPARAV.NLI ST.NL.mSEC .NEVAL.NFREO 
CSLTA , 0 EL‘ 4 IN'*TCL. 0 RATIG, ALPHA 
AFt ICl) , 3 F( ICI) ,CF( LCD ,CF( 10 l) 
AF, 3 F,C=, 0 P 


CO IDO I = 1 ,‘ILIST 
TE'-*o = ? APN ( I 1 


PARN( ! )=T = '^P*( l.*S!GN( I »*CELTA) 
IF(PAF N( I )-paRmaX( I ) ) L5» LC, IC 
PA=*i ( : ) = P ARMAXI I ) 

GC^O 25 

IC( piRvjl I )-PAR' 4 !M I ) ) 2 Cf 2 Cf 25 
DAR‘J { I ) = 3 ARmin( I ) 

CALL EVAL { =MN ,?ASN ) 

NEVAL=NEVAL ♦! 

I F( F'<N-F'.<P ) 3(-fACtAC 
pv-Js =m;j 

GC'^D ICC 


iC SIG'.M I )=-SIGN( I ) 

PAP" ( I )=TEHP* ( 1 . ♦ S IGM I ) *CELTA) 
I = ( -APNi I )-PAP«AM I ) ) tC, 5 :. 5 C 
5 D PA = ‘l ( I ) = P ARmaxI I ) 

GC’C 7 D 

f:C I=(?.F' 4 ( I )-PAPMlN( I ) ) 65 , 7 C, 7 C 
: 5 P ARN ( I ) = PAR“ IM I ) 

7 C CALL EVAL ( t ^ A 5 N ) 

N'iVAL='lEYAL + l 

1 F( =u‘i- = m 9 ) 3 C,FC,=C 

,J • CUP==M*| 

GC’C ICC 

9 C PARMI)=TE*iP 

id: cc'itin'Je 

: E’’J=N 
END 

5 i. 3 F'’L”!^'E -ancx :x,:y, = n:) 


■ L'-iEE^ Cr‘E=A’ 2 = 
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CCX -ICN S IGN' ( f : ) 50) » PARV IfS!( 50 ) 

CCMMON ICC} , ITYPEC 50) ,LIST(50) 

CCMMGM ,PCENTP ,.= START,FCEL ,GAIN0 

CC’^MCM NPiPAM ,NLI 5T f NLMSEC »NSVALrNFREQ 

CCMmcn C5LTA ..OELMIN ,TC L »OR A TI C , ALPHA 

complex AF{ ICl) ,BF( LCl ) ,CF{ LQl) ,0F(10l) 

CCMMON . AF,3F,CF,0F 
IY=IX*65529 
IF(IY) IC,20,K 
10 IY=IY+Z1A7A834A7+1 

20 RNO=IY 

PNOsP'lO^t. A 6 56612 2-9 

return 

EN-0 

SU3=GLri'jE =VAL( .3iS ) 

C SJFSO'JTINE E'.allaTES =IGLRE CF « = PI”. 

r 

c [men SI ON Pi=i:c) 

CG'^-'ON HGNI 50 ,PARMAX(50) f FA5MIN(50) 

CCMyGn PARAM(ICC) ,ITYPE{50) ,LIST(50) 

COMMON FPLVP ,FC£N^B ,rSTART,FCEL ,GA!NO 

CC^'^uN NPiPAy,NLlST ,NU**S=C ,NEVAL,NFPEC 

cgmmon celta ,oelmin,tcl,cratic,aloha 

CC'*PLEX AF( ici) ,5F( ICl) ,CF{ ICl) .DFIlOl) 

CCMmp.n AF,8F,CF,0F 

CC 5 I=1,MLIST 
IST( I ) 

5 PIRA-^ (L )=P AS ( I ) 

F“=0. 

EREi=FSTART 

:: :o i=i,-ifreo 

11= I 

CALL F ES=GN( GA IN ,F=EC ,I I ) 
t = -p=gain-ga:n3 

Fv»s CM >7 ’£mp«tcmp 

10 =F = ;= =REOi-=DEL 

FM=SO=Tf =M/N = 3=0 ) 

= E~ J.RN 
cN C 

SU3=3U7lNr GP AO H ( > ,C X , Y , N P ) 

T-I3 SJE=rUT!N= PLCTS F = ECLENCY RESPCi'SE. 

> f i c n 

LINE! 51)/l“t. ,lM./ 

SIGN) 5 r ) »PA P'<AX(5C) t^AP'^INMEO) 



o ri o 
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CC'^MGN PASA«( ICC)., ITYP = <50) ,LIST(50) 

CCM*«GN FPU'*P fFCEN'^P ,F STA5T ,FC=L ,GAINC 

NPAPA?* fNLIST.NLVSiC , N EVA L , NFP = 0 
CC.MMON DELTA ,OEL'* I N , TCL ,CS A TI C , ALPHA 

CC'^PLrX AF{ ICl ) ,3F( ICl) ,CF{ ICl) tCFdJl) 

CCMMCM AF,6F,CF,0F 

IWa 6 

tap IT=( U,666) 

666 FOR.MAT(lhl) 

YMAX»Y( 1) 

YMIN=YMAX 
CQ 10 1*2, MP 

IF ( Y( I ) .GT.YMAX ) YHA X=Yt I ) 

10 IF (Y( I ) .LT.Y^tn ) YMIN = Y(I) 

CI=Fa YMAX-YMIN 

I.= (OI‘==.N£.JJ 3CTC 2C 

nF IT = ( Ita, 300) Y^<A> 

300 =C5*< AT ( IX, • FOP all VALLES OF X, Y ECUALS dlPELZ.E) 

C-CTO 30 

20 *«5 iT = ( Ita, 100 ) 

1:0 FC = '^AT ( It— ///29X,c 1 ( • . • n 

IF(Y'^ax.LT.C.CP.Y'*I''I.GT.C) goto 5C 

I2=90.0«(-YMIN)/CIFF+1.5 

TrMPZ=L I.ME( IZ ) 

-L I.NEl IZ ) = V£PT 
50 DC cO 1=1, NP 

IY=90.:«(Y( I)-Y'^TN)/0IFr*1.5 
T?MP=L INE ( I Y) 

LINE( I Y )=STAP 
XX=X*( I-D-OX 

taP IT = ( Ita,20C) XX,Y( I ) ,LINF 
200 FCC'^A^ (IX, 2G12.5,^X,51A1) 

60 L !*! = ( lY )=TEHP 

/<? IT= ( IW, 400) 

40 0 =C=-^AT ( 2FX,ci( • . • ) ) 

JF (Y'" AX.G=.C.A\0.Y''!N.LE.C) LI.NE ( I Z) =Tr'-°Z 
3 0 5ET'J=N 

ENO 

SU3PCUTIME QLTPL7 

This SJ3S0UTIN5 PPI.V^S OLT FPECUENCY pespcsse. 

CCMVLN SIGN( fC) ,PAOFAX(5C) ,FAi:mI‘J(50) 

CC'^VCN PAP Au ( iCO) , ! TYPE ( 50 ) ,LI ST( 50) 

FPL'^o, >=C ENTS ,r START, FOE L,GAINC 
NPAPA'',^LIST,.^L'^SEC,^ = VAL,^FR = 0 
C:LTA.0EL'^IN,''CL,0=ATIC,AL=ha 
CO'-.^LEX A.= ( 101) ,3F{ LCD ,C=( 101) ,0^(101) 




<1 O n t** 
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OF POOR QU/ALlTY 


ccmm::n af,bf,cf,ci= 

I 6 

WP IT5I Iw, 600) 

600 FGP^'ATC////, l)«,LSFFP=gLcNCY PE SPC.NSE , // , 

> 19H FSEQUENCY GAIN,/, 

19H (GHZ) (Cei,//) 

R; EQaFSTART 
CO 20C I=1,NFREQ 
1 1* I 

CALL R£3PON(GA IN ,FRE0 ,II ) 

WRITE! IW. 610) FREC.G4IN 
610 FORMAT! 2F10 .3 ) 

2 30 F3=0=FR£0+f.j = l 
P ETURM 
EM 0 

SUEPOU'^INE OI£oly!FMP ,nIT) 

C this SJ5PGUTINE CISPLA>S CPTIMIZEO PASAVETE? 

/• 

V. 

Cl'^ ENSIGN 
CO‘^‘^GN 
COMmcn 
CCMMCN 
CC’’’*GN 
CG'^'^CM 

complex 

COMMON 
I w= 6 

W= IT = ! I,«, 600) M T.NEVAL ,F«D 
ibOO FOR*^ A"^ ( ///, 23H NC . CF ITERATION = .15,/, 

> 23H FLNCTICN5 EvALLA'^EC = .15,/, 

♦ 23H ERR0P=,E13.6) 

WP I'E ( !m , 61 :) 

610 PCR‘-'AT!////,21H OPTI'^IZED = ASAVE T = RS , // , 

> 23H SECT T\PE ? A S A '<E T= R $ , / / 


BUFFER! 1 C) 

SIGN! 50 ,RARV 4 X! 5 C) ,PARMIN! 50 ) 
P^RAM! ICC) , ITYPE! 50 ) .LIST! 5 0 ) 

= FLVP .FCENTR .FSTART ,FOEl ,GAING 
NPAPAV ,NLIST,\L.MS£C ,NEVAL, NFRSa 
C EL ^ A ,C ELf' IN .TCL.DP ATIC.AL PH A 
AF! ICD.EF! lCl),CF!i:i) ,CF! 131 ) 
AC,3F,CF,D = 


ENTRY POINT 


3 10 


entry GISP 

JLIST= 1 
IT=m?= 0 

CC AOC I = 1,NU'^SEC 
I = ( JL IST-NL 1ST ) 3CC 
J= I •^yf E( I ) 

GC'O (:lC,3iC,:lC,3 
N =3 


3 c : , A I : 


C ,32C ,320 ,320 ,330 ,3-0 


o n n 



GCT3 330 
330 N*5 


GOTO 35C 
340 N» I 

350 NMAX=ITcVP+N 
JFLAG«0 

C WAS THIS 5LE.M6NT VARIABLE 
355 IF(LIST(JLIST)-NMAX) 36C,36C,370 
360 JLIST»JLIST>L 
JFLAG= I 

IF(JL IST-nlIST) 355,355,37c 

WAS =LAG SET CURING ThE SE^aCH 

370 ITEwP-NMix 

IF(J=LAG) 4CC,4CC,3eC 
330 115*^?= I tEMP-aj 

CC 39 0 K=1,.N 
ITE'^P= ITEHP-l 
390 eUF = =:; (K ) = PAP AMI ITE'^F) 

WP I-^EI IW,640 ) I ,J ,(BLFF = P (L) ,L = l , N) 

640 FCaMATI I 4, 14, 5 = IC .4) 

4C0 ccntimue 
410 RETU3N 
ENC 

SU9R0LTINE SENSIT 
C 

r 

C This SUB = 3UTINE P£3 = CR**S sensitivity analysis. 

n 

r 

v» . 

CO'-MCS SIGN! fC) ,?AP ''AXI 5C) ,=AR'»IN(50) 

:G'-‘'^0N ?ARA«{ ICC) , I 7Y°E(50) ,LIST(50) 

cg'-*‘-':n fol«= ,fcentr ,f start, c^=l ,gaivj 

CC'^‘^ON N0iRAV,NLIST,NL*'SEC,NEVAL,N = = 5C 

CC'^'^ON CELTA.CEL'^IN ,TCL,DRATi:,AL3hA 

complex 4={ icn ,3F( ici) ,CF( LCD ,CF( 101 ) 

CCmmtnj A.= ,BF,CP,CF 

I W = 6 

s5 IT = ( I«,60: ) 

600 =2 5«AT(lhl,///,2S(ll-*),/, 
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I 

I 38H PRcOLENCY GAIN GAIN-P GAIN-**, 

I ♦ 24h S5NSIT-P SENSIT-f,/, 

' + 2X, 5H( GHZ ) f 6X ,4H| 08 J »A.M (09) ,6X ,A-H( 08) , 

, ♦ 5X,10H(P5P CENT) ,2X, ICHIPEP CEM),//) 

^ FREQ = rSTART 

t CO ^00 JalfNrPEC 

['■ JJ = J 

• TEHP=PARAM( I ) 

; CALL R5SPQN(Qej ,FREQ ,JJ) 

PAPAMI I )=1.02*TSr'P 

I CALL RESPGNIOBJP ,FPEC , JJ) 

r PARAM( I )=C.ga*TEHF 

I CALL a ESPON(OBJM ,FPEC, JJ ) 

SEMTYP* 5C.*(OSJ-C9 JP ) /C9 J 
! , S£NTY*^=50.«(Gej-C9JY ) /C9 J 

' '-R ITE( r«f 620) =REG,CS-,03JP ,C8J'^ ,'ENTYP,SEN'^Y'< 

i 62: FOR'*AT(Fi:.3,Fe.3,2FlC.3,2FL2.2) 

PARAM ( I ) = T£MD 

, FRE0= =REC+.=0EL 

AO: CCNTIN’JE 

I 50 J continue 

I OE'^UPM 

} END 
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APPEN13IX 3 APPROXIMATE EVALUATION OP DIODE TEST MOUNT 


EQUIVALENT CIRCUIT ELEMENTS 


Numarical values for the diode test mount equivalent circuit elements 
in Pig. 2.11(b) will be evaluated hers. Tho physical dimensions In 


Pig. 2.11(a) are as foUoKs: 


d = 2.03 mm 
h = 0.71 mm 


= 3*0^ mm 

D =7.00 mm 
o 


3.1 SvaluaMon of C^ 


The fringing capacitance, C^, Is given by Getslnger [103] as 


C^= 


r_2_ Jl.) 
^ b ’ b 


(3.1) 


Por the sxructuro in Pig. 2.11(a), 


D-d. 

0 


(3.2) 


D-d 

0 


(3.3) 


Substitution of numerical values into Sqs. 3.2 and 3.3 yields 


= 0.297 
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